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ABSTRACT 

Scattering  of  an  incident  plane,  or  cylindrical,  wave  by  a  staggered 
parallel-plate  waveguide  is  examined.  Of  particular  interest  are  the  modal 
coefficients  of  the  unattenuated  modes  in  the  waveguide  when  the  incident 
field  is  an  E-wave.  Several  past  approaches  based  on  ray  techniques 
are  critically  examined.  It  is  fcund  that  they  are  not  applicable  when 
the  stagger  is  small  or  when  the  incident  field  is  a  cylindrical  wave. 

For  the  latter  incidence,  we  introduce  here  an  aperture  field  method,  in 
which  the  modal  coefficients  are  calculated  through  the  Fourier  decomposition 
of  the  scattered  field  at  the  waveguide  aperture.  To  determine  this 
scattered  field,  we  apply  the  uniform  asymptotic  theory  for  the  primary 
contribution,  and  the  spectral  theory  for  the  first  interaction. 

Extensive  numerical  data  show  that  the  aperture  field  method  gives 
accurate  results  in  canonical  problems  whose  exact  solutions  are  known. 

As  an  application  to  the  design  of  the  constrained  waveguide  lens  antenna, 
our  method  enables  one  to  determine  the  field  distributed  over  the  lens 
in  a  more  exact  manner. 
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BASIC  CONVENTIONS 


1:  MKS  units  and  e-iwt  time  variation  are  used. 

2;  Y  -  (a2  _  k2)1/2  -  -i(k2  -  a2)1/2,  where  a  -  a  +  it  and 
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I .  INTRODUCTION 
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This  thesis  considers  the  problem  of  the  scattering  of  a  plane  or 
cylindrical  wave  by  the  staggered  parallel-plate  guide  in  Figure  1  when  the 
stagger  dimension  "l”  is  of  moderate  size  compared  to  the  wavelength.  In 
particular,  we  are  interested  in  finding  the  modal  coefficients  of  the 
unattenuated  modes  when  the  incident  field  is  an  E-wave.  The  problem 
goes  back  to  the  early  1940's  when  Kock  [18]  and  Ruze  [28]  realized  that  a 
waveguide  can  serve  as  a  refracting  medium  and  wide  angle  scanning  antennas 
can  be  designed  very  much  the  same  way  as  natural  dielectric  lenses.  While 
these  designs  proved  successful,  researchers  immediately  understood  that  for  a 
thorough,  quantiative  design  the  scattering  of  a  plane  wave  by  the  guide,  which 
constitutes  the  basic  element  of  the  lens,  had  to  be  well-understood.  Recent 
developments  in  satellite  communications  have  revived  interest  in  the  waveguide 
lens  [13],  and  make  the  problem  considered  here  all  the  more  relevant. 

Among  the  first  to- consider  the  scattering  of  a  plane  wave  by  a  guide 
have  been  Carlson  and  Heinz  [10],  Lengyel  [22]  and  Whitehead  [32].  However, 
they  studied  the  problem  in  which  the  parallel-plkte  guide  was  but  the 
single  element  of  an  infinite,  periodic  array  of  guides.  In  such  a  case, 
the  Wiener-Hopf  Method  and  other  classical  techniques  [32]  gave  an  exact 
solution.  These  solutions  applied  to  the  single,  isolated  guide,  however,  give 
poor  results  and  the  isolated  guide  has,  generally,  no  exact,  closed-form 
solution.  The  one  exception  to  this  statement  is  when  the  stagger  "J.  "  is 
zero.  Then,  the  Wiener-Hopf  Method  may  be  used  and  the  solution  obtained 
in  the  manner  outline  in  Chapter  II.  While  this  solution  is  useful  in  its 
own  right,  it  casts  little  light  on  how  to  proceed  when  is  not  zero  or 
very  small. 
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The  introduction  of  ray  techniques  to  the  study  of  scattering  by 
Keller  [17],  and  later  improvement  to  the  theory  of  Aklauwalia,  Lewis  and 
Boersma  [1]  provided  for  a  new  way  to  attack  the  problem.  Instead  of 
obtaining  the  full  solution,  rays  can  be  used  to  follow  each  scattering 
't  the  plates,  anc  from  each  scattering  the  modal  coefficients  computed 
by  Fourier  series  analysis.  This  approach  was  first  presented  by  Yee, 
Keller  and  Felsen  in  1968  [30];  a  critique  of  this  method  is  given  in 
Chapter  III.  In  Chapter  III,  we  show  that  the  YKF  method  presents  a 
^ i  ficult  singularity  problem  when  "z"  is  not  zero  and,  furthermore,  makes 
use  of  unwarranted  assumptions  about  the  ray  nature  of  the  scattered  fields 
in  computing  the  effects  of  multiple  scattering  between  the  plates.  A  much 
more  rigorous  use  of  ray  optics  for  the  unstaggered  guide  is  given  by 
Boersma  [2],  [5],  but  his  solution  is  not  easily  extended  to  the  staggered 
guide.  The  Modified  Diffraction  Coefficient,  given  by  Lee  [19],  [20],  is 
a  new  way  to  account  for  the  multiple  scattering  in  toto,  but  the  MDC, 
which  combines  functional  analysis  and  ray  optics,  can  be  properly  used 
only  when  "l"  is  large  -  usually  larger  than  one  wavelength.  Furthermore, 
all  the  previous  works  mentioned  above  examined  the  incident  plane  wave 
case  only. 

The  approach  to  the  problem  in  this  thesis  is  to  first  find  an  accurate 
approximation  of  the  total  field  on  the  guide  aperture,  and  then  by  Fourier 
decomposition  compute  the  modal  coefficients.  Because  of  the  difficulty  in 
formulating  the  multiple  scattered  fields  when  "i"  is  not  zero  or  large 
we  include  in  the  aperture  field,  for  the  incident  plane  wave,  only  the 
first  interaction,  or  scattering,  between  the  plates  of  the  guide.  Previous 
plane-wave  investigations  of  the  scattered  field  that  included  the  first 
interaction  were  conducted  by  Jones  [16]  and  Boersma  [4],  but  their  results 


respectively  were  in  terms  of  very  el&borate  doubly  infinite  integrals  and 
series.  In  this  thesis  the  primary  scattering  field  is  given  by  the  canonical 
or  UAT  solution  to  the  scattering  by  a  half-plane,  and  for  the  plane  wave 
the  first  interaction  is  found  by  the  Spectral  Theory  of  Diffraction,  recently 
introduced  by  Mittra  and  Rahmat-Samii  [26].  This  formulation  of  the  aperture 
field,  given  in  Chapter  IV,  is  relatively  easy  to  compute  and  is  equally  valid 
for  all  stagger  size  and  plate  separation  "a."  Using  this  aperture  field, 
extensive  data  for  the  modal  coefficient  of  the  important  TE^q  mode  are 
presented  in  Chapter  IV,  and  an  application  to  the  study  of  waveguide  lens 
antennas  is  given.  Finally,  conclusions  are  presented  in  Chapter  VI. 

Finally,  we  note  several  conventions  used  in  this  thesis.  (1)  The 
time  factor  is  e  and  is  suppressed.  (2)  All  angles  are  between  -ir 
and  fir,  positive  counterclockwise,  unless  otherwise  noted.  (3)  The  E-wave 
problem  is  examined  in  detail,  the  extension  to  the  H-wave  problem  is  only 
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II.  SUMMARY  OF  THE  WIENER-HOPF  METHOD,  THE  SPECTRAL 
THEORY  OF  DIFFRACTION  (STD)  AND  THE  UNIFORM 
ASYMPTOTIC  THEORY  OF  DIFFRACTION  (UAT) 


For  convenience  of  reference  and  as  an  introduction  to  the  notation 
used  in  this  work,  the  Wiener-Hopf  Method,  the  Spectral  Theory  of 
Diffraction,  and  the  Uniform  Asymptotic  Theory  of  Diffraction  are  briefly 
summarized  below. 

A.  The  Wiener-Hopf  Method 

Before  we  come  to  the  problem  of  scattering  by  a  half-plane,  a  quick 
discussion  of  some  elements  of  the  Fourier  transform  complex  number  theory 
relevant  to  the  Wiener-Hopf  Method  is  presented. 

A. 1.  Mathematical  preliminaries 

First,  if  <{i(x,z)  and  <t(x,a)  form  a  Fourier  transform  pair,  with 
a  -  a  +  it  and  if  <j>(x,z)  has  an  asymptotic  behavior  given  by  <t>(x,z)  '  AeT  z, 
as  z  -►  +»,  then  the  half-range  Fourier  transform,  <f+(x,a),  defined  by 


4>+(x,a) 


1 


<(>(x,z)eictZ  dz  , 


(II. 1) 


will  be  analytic  in  the  region  t  >  t  .  This  statement  is  true  for  any 
<P(x,z)  that  is  "well-behaved"  and  for  any  kernel  that  is  continuous  in 
a  and  z  and  analytic  in  z.  In  II. 1,  it  is  eictz.  A  more  rigorous  formulation 
of  the  sufficient  conditions  involved  in  this  discussion  can  be  found  in 
[25]  or  [29].  We  assume  throughout  that  $(x,z)  is  such  that  the  conditions  in 
[25]  or  [29]  are  indeed  satisfied  and  the  final  results  may  be  applied  in  a 
straightforward  manner. 

+ 

Likewise,  if  4>(x,z)  ~  Be1  Z,  z  -*■  -«,  then  $_(x,a),  which  is  given  by 

$_(x,a)  - [  <Kx,z)eiaZ  dz  ,  (II. 2) 

y2tt  ' 


•A, .  . 


I 
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is  analytic  in  the  region  x  <  x+.  Note  that  the  subscripts  "+"  and  can 
now  be  used  to  indicate  that  the  function  involved  is  analytic  in  some 
"plus"  or  "minus,"  that  is,  upper  or  lower  region  in  the  a  plane. 

We  now  define  two  functions  4> (x, z)  and  ^  ^  (x,  z)  as 


4>(+^(x,  z)  - 

and 

<*>(_^(x,  2)  - 


<$>(x»  z)  z  >  0 

0  z  <  0 

r  0  2  >  0 

,$(x.  z)  2  <  0 


(II. 3a) 


(II. 3b) 


Then,  from  the  analytic  properties  of  and  $+  in  their  respective  regions, 
it  is  easily  shown  that  <J+  and  $(  ),  <f  are  Fourier  transform  pairs. 

Also  essential  to  the  Wiener-Hopf  method  IS  the  factorization  and 
decomposition  of  complex-valued  functions.  By  factorization,  we  mean 
expressing  a  function  G(a)  as  the  product  of  "plus"  and  "minus"  functions, 
that  is,  G(a)  -  G+(a)G  (a),  where  G+(a)  is  regular  for  x  >  x  and  G  (a) 
is  regular  for  x  <  x+,  x  ,  x+  some  numbers.  For  a  well-posed  problem,  we 
will  have  x+  >  x  so  that  a  region  exists  in  which  G+(a)  and  G  (a)  are  both 
regular.  In  contrast,  decomposition  refers  to  expressing  a  function  S(a)  as 
a  sum  of  "plus"  and  "minus"  functions  such  that  S(a)  -  S+(a)  +  S  (a),  where, 
as  before,  we  expect  to  find  some  region,  x  <  x  <  x+,  in  which  S+(a)  and 
S_(a)  are  both  analytic. 

A  useful  formula  that  is  later  used  for  factoring  an  entire,  even 
function  G(a)  is 


G(a)  -  G(0)  "  (1  -  -2-)e-ia/an  S 

i-1  lYn  i-1 


(1  + 


a  .  +ia/an 


iY 


-)e 


(II. 4) 


n 


It  is 


where  +iy  are  the  simple  zeroes  of  G(o)  and  y  -*■  an  +  b,  n  °>. 

-  n  n 

obvious  that  the  first  product  is  analytic  in  the  upper-half-plane  given 
by  :  >  t",  where  x”  is  determined  by  the  location  of  Yq  in  the  a  plane. 

For  decomposition  of  a  function  S(ot)  analytic  in  some  strip 
x  <  t  <  x+,  and  which  has  the  property  that  |s(a)|  -*•  0  uniformly  as 
1  o |  -*•  «  within  that  strip,  we  will  have  occasion  to  use  the  following  formulas. 

“H-ic 

(ffil  de  i  <  c  <  t+  (II. 5a) 

-®+ic 


S+(o) 


2iri 


and 


®+id 


S  (a) 


JL 

“  2iri 


S(B) 
(B  -  a) 


dB  t  <  d  <  t 


(II. 5b) 


-®+id 


That  is,  within  the  strip,  using  Equations  (II. 5a)  and  (II. 5b),  we  can  write 
S(a)  ■  S+(a)  +  S_(a). 

The  above  equations  for  decomposition  are  also  useful  for  factorization. 
If  G(a)  is  analytic  and  nonzero  in  some  strip  a  <  t  <  b,  then  by  taking 
In  G+(a)G  (a)  we  can  use  (II. 5a)  and  (II. 5b)  to  find 


G+(a)  »  exp 


- 

®+ic 

- 

1 

In  G(cO 
(B  -  a)  dS 

2iri 

_ 

-oo+ic 

X  <  C  <  X 


(II. 6a) 


and 


G  (a)  ”  exp 


®+id 


1 

2iri 


-oo+id 


x  <  d  <  x 


(II. 6b) 


Now,  however,  we  require  that  within  the  strip  |G(a)|  -*■  1  uniformly  as 
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j a j  -*■«*.  Notice  that,  because  of  the  use  of  the  exponential,  G+(ct)  and 
G  (a)  are  also  nonzero  in  their  regions  of  regularity.  Other  factorization 
and  decomposition  formulas  can  be  found  in  [24]  or  [25].  Having  finished 
these  mathematical  preliminaries,  we  now  come  to  the  problem  that  illustrates 
the  three  methods  in  this  chapter  -  scattering  by  a  half  plane. 

A. 2.  Scattering  by  a  half -plane 

Consider  the  problem  in  Figure  2.  We  wish  to  find  the  scattered 
field  when  the  incident  field  is 


^(x  z)  -  eik(xsin£i+zcosft)  m  eikpcos(G-6) 


(11.7) 


where  <b *  can  be  for  an  E  wave,  and  for  an  H  wave.  Without  loss  of 

y  y 

generality,  we  place  the  restriction  it  >  Q  >  0.  There  are  now  two  ways  to 

start  formulating  the  Wiener-Hopf  problem:  (a)  use  the  Green's  function 

method  which  then  leads  to  an  integral  equation,  or  (b)  use  the  differential 

equation.  Although  the  Green's  function  method  is  more  general,  the  second 

approach,  called  Jones'  Method  [15],  is  simpler  and  clearly  shows  the 

requirements  for  the  Wiener-Hopf  solution. 

We  now  examine  the  E-wave  problem  in  some  detail  and  later  state  the 

t  s  is 

equivalent  results  for  the  H-wave.  First  let  $  *  $  +  $  ,  <f>  being  the 

scattered  field.  Then  $  satisfies  the  source-free  Helmhotz  equation 


32$s  +  d2<pS  +  kZ<t>3  -  0 
x  z 


(II. 8a) 


with  the  boundary  condition 


<(>S(0,z)  »  -$*(0,2)  -  »  <  z  <  0 


(II. 8b) 


'jgjm  - 


The 


A  small  but  finite  loss  is  introduced  by  letting  k  *  k^  +  ^2* 

Wiener-Hopf  solution  begins  by  transforming  (II. 8a)  and  (II. 8b)  into  the 
spectral  or  transform  domain,  by  multiplying  both  sides  of  (II. 8a)  and 
(II. 8b)  by  eiaz//2jT  and  integrating  over  all  z.  If  this  operation  is  performed 
Equations  (II. 8a)  and  (II. 8b)  then  become 

32$S(x,  a)  -  y2<J’S(x,  a)  -  0  (II. 9a) 

and 

$s(0,  a)  »  -0*  (0,  a)  t  <  -k2  cos  ft  ,  (II. 9b) 

where  the  subscript  "minus"  appears  because  (II. 8b)  is  true  only  for 

2  2  1/2 

-“>  <  z  <  0.  The  function  y  **  (a  -  k  )  is  defined  by  its  value  on  the 
upper  or  proper  sheet,  as  shown  in  Figure  3.  These  branch  cuts  insure 
that  the  radiation  condition  is  satisfied  for  all  a  on  the  real-number  line 
and  that  Re(y)  >  0  on  the  entire  sheet.  Note  that  y  is  analytic  in  the 
region  k2  >  r  >  -k2* 

The  solution  to  (II. 9b)  is 

$s(x,  a)  *  A(a)e“Y^  ,  (II. 10) 

so  that  at  x  *  0  we  have 

-$*(0,  a)  +  $®(0,a)  =*  A(a)  (II.  11) 

where,  since  satisfies  the  radiation  condition,  is  analytic  for 
t  >  ^2*  Defining  J_(a)  as  being  proportional  to  the  transform  of  the 
total  induced  current  on  the  plate,  we  get 

J_(c)  •  «)  -  V-*0"’  t‘)1  "  ~',kM  ■  (11.12) 


*  y, 


Vr  * 


r 
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Combining  (11.12)  with  (II. 11),  we  obtain  the  Wiener-Hopf  equation 

$®(0,  a)  +  J_(a)  -  $*(0,  a)  -  0,  -k£  cos  SI  >  x  >  (11.13) 

(see  Figure  4),  which  is  in  the  standard  Wiener-Hopf  form 

R(a)<J>^(a)  +  G(a)J  (a)  +  T(a)  =>  0  t+  >  t  >  t  ,  (11.14) 

«  g 

with  R,  G,  T  given  and  analytic  in  t  >  r  >  t  and  the  unknowns  $+,  J 
analytic  in  :  >  t  ,  t  <  t+,  respectively.  Although  Equations  (11.13)  and 
(11.14)  are  algebraic  equations  in  two  unknowns,  by  using  the  Wiener-Hopf 
Method,  a  solution  is  still  possible. 

The  first  step  is  to  factor  G(a)  *  (1/y)  into  its  "plus"  and  "minus" 
parts.  For  this  particular  case,  it  is  obvious  that 

G+(a)  -  (a  +  k)'1/2  ,  (II. 15a) 

and 

G_(a)  •  (a  -  k)~1/2  .  (II. 15b) 

Dividing  both  sides  by  G+(a) ,  it  follows  that 

4>®(0,  a)  (a  +  k)1/2  +  (a  -  k)1/2  J_(a)  -  (a  +  k)1/2  $*(0,  a)  =  0.  (11.16) 

1/2  i 

The  next  step  is  to  decompose  (a  +  k)  $  .  Again,  the  decomposition 
can  be  done  by  inspection,  and  with  $^(0,a)  given  by 

<P^(0,  a)  *  -i/[/27f(a  +  k  cos  55)],  t  <  cos  55  ,  (11.17) 

one  can  write 


A 

V 


'  A 


t*  :*<&■  j* 


/ 


i(g  +  k) 


1/2 


/lit  (a  +  k  cos  ft) 


_  If (a  +  k) 1/2  -  (k  -  k  cos  ft)1/2 
/in  {a.  +  k  cos  ft) 
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1/2 

1  (k  -  k  cos  ft) 

^2^  (o  +  k  cos  ft) 

The  first  term  on  the  right-hand  side  is  S+(n)  which,  because  of  the 


(11.18) 


cancellation  of  the  singularity  at  a=  -kcos  ft,  is  analytic  for 
t  >  -k2.  The  second  term  on  the  right  is  S  (o) ,  and  it  is  obviously 
analytic  for  t  <  -k2  cos  ft.  Combining  (II. 15a)  thru  (11.18),  we  have 

(«  +  k)1/2  *?(o,  g)  -  *•! fc. ..t-KL17 2 . A  -  -cos  Q) 1/2 

/2?(a  +  k  cos  ft) 


J  (g) 


(g  -  k) 


1/2 


1/2 

l(k  -  k  cos  ft) 
/2n(a  +  k  cos  ft) 


(11.19) 


for  -k2  cos  ft  >  x  >  -k2-  The  right  side  of  (11.19)  is  analytic  in  the 
upper-half  plane,  t  >  -k2>  the  left  side  analytic  in  the  lower-half  plane, 
t  <  -k2  cos  ft.  Since  there  is  a  region  where  both  sides  of  (11.19)  are 

analytic  and  equal  to  each  other,  it  is  possible  to  analytically  continue 
both  sides  of  (11.19)  into  the  entire  complex  plane,  and  obtain  the  equations, 

(o  +  k)1/2  *®(0,  g)  -  i^ot.  +-k2-  / — c-os  a)-  ' ■  =  p(o)  (II. 20a) 

/27(g  +  k  cos  ft) 


and 

J  («)  0a/2 

.  _ = _ i(k  -  k  cos  ft)  „  p,  . 

i/2  . —  » 

(g  -  k)  /2rT(g  +  k  cos  ft) 

where  P(g)  is  an  entire  function.  We  see  that  the  Wiener-Hopf  Method 


(II. 20b) 
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leads  to  two  equations  in  which  the  two  unknowns  may  be  solved  for 
within  some  entire  function,  P(a). 

To  completely  solve  (II. 20a)  or  (II. 20b),  another  condition  must  be 
applied  to  find  P(a).  In  electromagnetics,  this  condition  is  the 
"edge  condition,"  [23],  which  states  that  no  charge  may  accumulate  at 
the  edge.  This,  in  turn,  implies  that  $“(p,  9)  *  0(p  ) ,  p  0,  1>6>0, 
and,  therefore,  $®(0,  a)  *  0(a  ^  ^)  as  (a|  -*■  «.  Inspecting  the  asymptotic 
behavior  of  (II.20a)  and  using  Liouville's  theorem  P(a)  must  be  zero. 

Then,  from  (II. 20a)  and  (11.12) 


J  (a) 


i(k  -  k  cos  0)1/2(o-  k)1/2 


/Hr  (a  +  k  cos  0) 


*yA(a) 


(11.21) 


Finally,  solving  for  A (a)  and  inverse  Fourier  transforming  $S(x,  a)  we  obtain 


the  desired  result 


(x,  z)  - 


CO  p 

[  2(k  -  k  CO; 

i  L  + 


£°-k-^1/-2^-  da.  (11.22) 

+  k  cos  Q)  ,  2y 


s  s 

For  an  H-wave,  $  ■  Hv>  the  boundary  condition  in  the  transform  domain 
is  3x$^(0,  a)  »  -8x$*(0,  »  9S(x,  a)  -  sgn(x)A(a)  e  Y|X^,  and  the  current, 
which  is  polarized  in  the  z  direction,  is 


J_(a)  =*  -|{b®(0+,  a)  -  $f(Cf ,  a) }  ■  A(a)  . 


(11.23) 


If  the  same  considerations  as  for  the  E-wave  are  followed,  the  scattered 


field  for  the  H  wave  is 


s  m  3gn .(xj.  7  2 (k  4-  k  cos  fl)1/2(a  +  k)1/2  e~Y'X  -iaz 

’  2ir  J  (a  +  k  cos  SI)  2y 

—CO 


dz. (11.24) 


r 
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In  summary,  once  an  equation  of  the  type  (11.14)  is  derived,  the 
Wiener-Hopf  method  allows  one  to  solve  it  for  its  two  unknowns  within 
some  entire  function  P(a).  The  essential  features  of  the  Wiener-Hopf 
Method  are  then  factorization  and  decomposition,  the  use  of  analytic 
continuation,  and  the  necessity  of  imposing  another  physical  constraint 
to  obtain  a  full  solution. 

B.  Spectral  Theory  of  Diffraction  (STD) 

In  the  course  of  investigating  the  coupling  into  a  waveguide,  we 
need  to  consider  the  scattering  by  a  half -plane  of  rather  complicated 
incident  fields.  While  the  problem  of  a  plane-wave  diffraction  by  a 
half  plane  has  been  extensively  analyzed  since  Sommerf eld’s  solution  in 
1896,  the  diffraction  of  a  more  general  incident  field  has  only  recently 
been  examined.  The  first  treatment  of  an  incident  field  other  than  a  plane 
wave  was  by  Carslaw  in  1899  who  considered  the  very  special  case  of 
an  incident  field  due  to  an  isotropic  line  source.  More  recent  work 
with  the  isotropic  line  field  can  be  found  in  Born  and  Wolf  [8],  and  in 
the  book  by  Clemmow  [11],  who  examined  the  problem  through  the  use  of  a 
plane-wave  spectrum  representation.  A  high  frequency  or  asymptotic 
solution  to  this  problem  and  of  a  more  general  line  source  field  was 
studied  by  Boersma  and  Lee  employing  the  U.A.T.  formalism  [6],  [7]. 

The  first  rigorous  and  extensive  study  of  the  case  of  an  arbitrary 
incident  field  scattered  by  a  half-plane  was  done  by  Rahmat-Samii  and 
Mittra  [26],  and  their  approach,  labeled  the  Spectral  Theory  of  Diffraction, 
is  outlined  below. 


We  start  by  examining  in  more  detail  the  expression  for  the  scattered 
field  of  an  incident  plane  wave  as  given  by  Equations  (11.22)  and  (11.24). 
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Both  E  and  H  waves  can  be  conveniently  considered  by  letting 


'-1 

E 

wave 

+1 

H 

wave 

(11.25) 


The  "spectral  diffraction  coefficient”  is  defined  as 
o)  ,  2(k 

2  (a  +  k  cos  SI) 

where  the  square  roots  are  defined  by  taking  the  cut  just  below  the 
negative  real  axis  in  the  a  plane.  The  scattered  field  can  therefore 
be  now  written  as 


(11.26) 


'  1  ■ 

i  r®  -ylxl 

E  wave 

4>S  (x,  2)  * 

sgn(x) 

2tt  j  X(V  0)  2y  e  da’  f0r 

mOO 

H  wave 

(11.27) 


Note  that  (e~Y^x'/2y)  is  the  Fourier  transform,  with  respect  to  2,  of 
i  /  o  9 

iHQ(k/x  +  2  )/4,  which  is  the  Green's  function  for  a  line  source 
radiating  in  free  space. 

Introducing  the  change  of  variables  x  *  p  sin  0,  z  ■  p  cos  9 
a  »  -k  cos  4>  and  y  *  -ik  sin  ip  into  (11.27),  one  arrives  at 


♦ (P,  9) 


’  1  ' 

sgn(0); 

4tt 

r 

x(n,  *)eikp  cos  (*  -  |9l)  d* 


(11.28) 


In  the  preceding  equation,  ip  is  the  complex  angle  defined  on  the  path 
T  shown  in  Figure  5,  and 


X(fl,  <P)  -  X(k  cos  fl,  -  k  cos  ^)  -  xi(J7,  <P)  +  exr(^»  P) 


(11.29) 


a 


where 


X 


i,r 


£  CSC 


'n  ±  ip' 

2 


(11.30) 


Clearly  and  Xr  are  infinite  at  ip  ■  0  and  at  <l>  ■  -Q,  respectively. 

These  two  values  of  ip  correspond  to  the  incident  and  reflected  shadow 
boundaries  appearing  in  the  GTD  technique.  As  a  matter  of  fact,  the 
spectral  diffraction  coefficient  is  precisely  the  angular  part  of 
Keller's  diffraction  coefficient,  when  tp  is  replaced  by  the  observation 
angle  0.  Although  x  tends  to  infinity  at  the  shadow  boundaries,  the 
scattered  field,  as  evident  in  (11.27),  does  not. 

From  (11.27)  and  the  fact  that  (e_Y^/2y)  is  the  Fourier  transform  of 
the  Green's  function  (iHg(kp)/4),  it  follows  that  x(^»  0  is  proportional 
to  the  Fourier  transform  of  the  induced  current  on  the  plate.  We  may, 
therefore,  define  xP°(n>  <10  as  the  transform  of  the  physical  optics 
part  of  the  current  and  write 


XP°(fl>  <i>)  *  XP°(n»  0  +  exP°(ft,  0 


(11.31) 


where 


po 

Xi,i 


*  +  ctn 


n  +  it; 


Then  x^defined  by 


Xf<n>  0  m  x(n,  ip )  -  xpo(n»  0 


(11.32) 


(11.33) 


is  bounded  everywhere  and  could  be  called  the  spectral  fringe  diffraction 
coefficient.  It  is  used  in  scalar  aperture  diffraction. 


Consider  again  the  scattered  field  in  (11.27).  This  spectral  integral 


can  be  exactly  expressed  in  terms  of  the  Fresnel  integral,  viz., 

*s(p,  0)  -  _eikPcos(n-e)F(_ci)  +  eeikpcos(^0)F(cr) 


(11.34) 


(11.35) 


(11.36) 


So  far  in  this  development,  we  have  tacitly  assumed  that  n  is  real. 
However,  using  the  argument  of  analytic  continuation  one  can  show  that 
if  the  angle  of  incidence  is  a  complex  number  w.  Equation  (11.34)  is  still 
a  valid  solution  but  with  fl  replaced  by  the  complex  u  in  (11.34)  and  (11.36). 

With  the  solution  for  a  plane  wave  of  complex  incident  angle  available, 
the  case  of  an  arbitrary  incident  field  on  a  half  plane  as  shown  in 
Figure  6  can  now  be  examined.  We  assume  the  incident  field  can  be 
placed  in  a  spectral  representation  as 

*  ,  f  ikp.  cos  (oj—  !  0_  | ) 

♦  (Pl,  0j>  -  P(u)e  dw  (11.37) 

r* 


where  p^  «  0  denotes  the  phase  center  of  the  incident  field  and  the  path 
T'  is  the  same  as  the  path  T  but  where  I"  accounts  for  any  singularities 
in  P(oj).  Expressing  now  in  terms  of  the  coordinates  attached  to  the 


♦A(P,  9) 


i-  f 

4  if  J 


[P(“)e 


ikdcos 


(u>-ft) .  ikp 
Je 


cos(w-9) 


doj 


(11.38) 


when  p  sin  $  >  -d  sin  ft.  The  incident  field  is  then  a  spectrum  of  plane 
waves  incident  at  a  complex  angle  ai  with  a  weight,  or  spectrum  density, 
given  by  the  term  in  the  brackets.  Since  (11.34)  is  the  solution  for  a 
plane  wave  incidence,  even  with  a  complex  angle  of  incidence,  by  super¬ 
position  the  scattered  field  due  to  an  arbitrary  incident  field  is 


4>s(p ,  0) 


_L_  (  [P(a))eikdcos(“-fl>]{eikPcos(“-9>F(_ci) 

4fri  J 

r 


-  £e±kpc0s(ari-e>F(5r)}da1  (11.39) 

where  F  is  the  Fresnel  function  of  a  complex  argument,  discussed  in 
1  IT 

Appendix  A,  and  5  *  are  complex  valued  and  again  given  by  (11.36)  with 
ft  replaced  by  oj.  We  make  use  of  (11.39)  in  Chapter  IV. 

C.  Uniform  Asymptotic  Theory  of  Diffraction  (UAT) 

Except  for  a  few  canonical  geometries  commonly  found  in  textbooks, 
closed-form  solutions  for  Maxwell's  equations  are  not  known.  Several 
asymptotic  or  high  frequency  techniques  have  therefore  been  devised 
that  enable  one  to  study  the  scattering  process  locally,  thereby  making 
good  use  of  the  canonical  solutions.  One  of  these  high  frequency  theories, 
the  Uniform  Asymptotic  Theory  of  Diffraction  (UAT) ,  proposed  by 
Ahlawalia,  Lewis  and  Boersma  [  1  ] ,  is  based  on  the  form  of  the  scattered 
field  in  Equation  (11.34).  The  UAT  is  one,  and  the  most  successful, 
of  several  high  frequency  or  asymptotic  theories  that  attempt  to 
overcome  the  singularity  of  the  diffracted  field  as  formulated  by  the  GTD  [7], 
For  the  canonical  problem  of  an  incident  plane  wave  on  a  half -plane 
(Figure  7),  all  asymptotic  solutions  reduce  to  the  one  in  (11.34). 

But  for  a  more  complex  situation,  such  as  a  line  source  field  incident 


Vr 


F 
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on  a  half -plane,  it  has  been  shown  [7]  that  the  UAT,  which  provides  a 
systematic  way  for  obtaining  the  higher-order  terms  in  the  asymptotic 
expansion,  gives  these  terms  accurately,  whereas  the  other  methods  do  not. 
Furthermore,  the  UAT  has  been  successfully  applied  to  problems  with  a 
patterned  incident  field  and  for  more  complicated  geometries  of  the 
scatterer  [12],  [26], 

The  Uniform  Asymptotic  Theory  requires  that  the  incident  field  be 
a  ray  field,  that  is,  $*(r)  can  be  expressed  as 


♦hr)  ■  elks(r>[J(?)  +  l  (lk>-*  (?>] 

n-1 


(11.40) 


where  s(r)  is  the  eikonal  and  the  A’s  are  amplitude  vectors.  Neglecting  in 
(11.40)  the  k  1  and  higher  terms  and  requiring  that  iji1  satisfy  Helmholtz's 


equation. 


(V2  +  k2)*1  -  0  , 


(11.41) 


we  have 


k2(l  -  (Vs)2)A  +  ik(AV2s  +  2Vs  •  VA)  +  VA  »  0 


From  Equation  (11.42),  the  following  conditions  are  derived: 


(11.42) 


IVsT  -  1 


2Vs  •  VA  +  (Vs) A  -  0 


and  from  Gauss'  Law 


(II. 43a) 

(II. 43b) 


1 


Vs  •  X  *  0 


(II.43c) 
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Also  from  (11.42),  we  must  have 


«  0 


(11.44) 


or  very  close  to  it.  This  last  condition  obviously  does  not  hold  at  a 
caustic  or  at  other  points  where  the  field  variation  is  great. 

If  the  incident  field  is  indeed  a  ray  field  and  satisfies  the  previous 
conditions,  then  the  solution  to  the  scattering  problem  shown  in  Figure  7 
is  given  by  UAT  as 


d«t(r)  -  <{iG  +  4id  +  0(k"3/2) 


where  4>a (r)  is  called  the  diffracted  field  and  is  given  by 


<j>d(r)  -  <fi(r  -  0)D(ft,  0)  g(kr) 


(11.45) 


(11.46) 


where 


g(kr) 


ik(r+-ir/4) 


(11.47) 


4  cos  (-|ft)  sin  C^Q) 

and  D  -  - - - 3 -  for  H  waves 

cos  ft  -  cos  9 


4  sin  (-|ft)  cos  Cj0) 
cos  ft  -  cos  0 


for  E  waves 


(11.48) 


The  angles  are  illustrated  in  Figure  7.  Note  that  at  the  lines  0  *  +  ft 
the  diffraction  coefficient  D  becomes  infinite  and  the  diffracted  field 
diverges.  These  lines  are  called  the  incident  and  reflected  shadow 
boundaries . 


- taster 
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The  first  term  in  (11.45),  <$>  (r),  which  represents  a  modification  of  the 
geometrical  optics  field,  is  given  by 


$G(?)  -  [FCC1)  -  FCS1))*1-^)  +  [F(£r)  -  F(£r)]*r(r) 


(11.49) 


For  large  £,  the  Fresnel  function  has  an  asymptotic  expansion  given  by 


F(£)  -  Q(-£)  +  F(£)  +  0(£"3)  |£j  -  « 


(11.50) 


where  0(x)  is  the  unit  step  function,  <£r(r)  is  the  field  radiated  in  free 


space  by  the  image  source  of  $  (r),  and 

i«>  .  -L_  et«2«/4)  . 

2£/tT 

The  term  £*,r  is  called  the  detour  parameter  of  the  incident  (reflected) 
field  and  is  a  measure  of  the  phase  difference  between  the  incident 
(reflected)  and  diffracted  fields.  Specifically  (refer  to  Figure  8), 


(11.51) 


£i*r(?)  -  e1’*  l/kU  +  b-c1’")! 

•{  y» 

where  £  ’  (r)  is  the  "shadow  indicator"  and  is  defined  as 


(11.52) 


e  (r) 


i,r 


+1  if  r  is  in  the  geometrical  shadow  of  <j>  * 
-1  if  r  is  in  the  geometrical  lit  region  of 


(11.53) 


Observe  that: 

i  r 

(a)  If  £  *  is  large,  that  is,  if  we  are  away  from  the  shadow  boundaries, 
Q  2 

<t>  reduces  to  <f>°,  the  geometrical  optics  field.  The  total  field 
is  then  given  by 


f>C(r)  -  <f>8(r)  +  $d(r)  . 


(11.54) 


r' 

•»**- 


j. 


'■ 
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This  form  of  Che  solution  is  known  as  the  GTD  (Geometrical 
Theory  of  Diffraction)  field. 

(b)  If  we  are  on  the  shadow  boundaries,  5*’r  becomes  zero  and 

F(5  ’  )  diverges.  However,  recall  that  on  the  shadow  boundaries 
the  diffraction  coefficient  diverges  as  well.  It  is  easy  to 
show  that  the  singularities  of  F(Ci,r)  and  D(D,9)  cancel, 
making  the  total  field  well-behaved.  It  is  in  this  sense  that 
the  UAT  is  "uniform." 

G  st 

To  find  the  region  in  space  where  $  approaches  ,  we  seek  to  find 
the  locus  of  points  such  that 

Ui,r|  *  M  (11.55) 

_3 

where  M  is  some  constant  large  enough  so  that  the  0(5  )  term  in 

Equation  (11.50)  can  be  ignored.  The  locus  of  points  satisfying  (11.55) 
is  obviously  a  hyperbola  with  one  focus  at  the  source  (image)  and  the 
other  at  the  edge.  The  asymptotes  of  the  hyperbola  form  an  angle  26, 
where  (see  Figure  9) 

6-4  sin-1/(M/  2ka)  .  (11.56) 

G  8 

For  M  -  2,  the  difference  between  $  and  is  approximately  2%.  For  M  -  1, 
it  is  approximately  10Z.  Note  that  for  small  M,  that  is,  in  the  "transition 
region,"  the  field  in  (11.45)  is  not  a  ray  field. 

Finally,  we  note  in  connection  with  the  UAT  that  if  the  line  source 
in  Figure  7  gives  rise  to  a  patterned  field,  or  if  higher-order  terms 
are  desired,  (11.40)  through  (11.44)  provide  the  machinery  for  doing  so. 
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III.  SOLUTIONS  OF  CANONICAL  PROBLEMS 

In  order  to  gain  insight  into  our  study  of  coupling  into  a  parallel 
plate  guide  with  an  arbitrary  stagger  "£"  (Figure  1) ,  the  case  of  the 
unstaggered  waveguide  (£  ■  0)  is  examined  in  some  detail.  This  canonical 
problem,  which  has  an  exact  closed-form  solution,  also  provides  the 
necessary  test  for  evaluation  of  the  different  staggered  waveguide 
solutions  presented  in  Chapter  IV. 

Because  of  the  symmetry  of  the  unstaggered  guide,  we  can  introduce 
an  electric  or  magnetic  wall  halfway  through  the  guide  and  study  the 
corresponding  cases  in  Figures  10a  and  10b.  The  solutions  to  the 
scattering  by  the  two  guides  in  Figures  10a  and  10b  now  correspond 
to  the  symmetrical  and  asymmetrical  components  of  the  total  scattered 
field.  We,  therefore,  speak  interchangeably  of  the  canonical 
case  as  being  a  guide  of  zero  stagger  or  one  where  the  stagger  l  is  infinite, 
the  latter  case  sometimes  allowing  an  easier  interpretation  of  the  results 
and  a  more  natural  extension  to  the  case  where  is  finite. 

There  are,  essentially,  two  approaches  to  solving  the  staggered,  or 
equivalently,  the  infinite  staggered  case.  The  first  is  to  apply  some 
function  theoretic  technique  and  derive  the  modal  coefficients  in  the  guide 
with  the  interactions  between  the  plates  accounted  for  in  toto.  This  approach, 
using  the  Wiener-Hopf  Method,  is  illustrated  in  Sections  A  and  B  wherein 
the  "ray-to-mode"  conversion  formulas  and  the  "modx^ied  diffraction 
coefficient"  are  presented.  In  Section  C,  we  look  at  a  ray  technique 
the  YKF  Method  [30],  which  tries  to  track  the  individual  interactions 
between  the  plates,  and  present  this  method's  "ray-to-mode"  conversion 


1 

I 
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Figure  II.  1)1  f f raci  ion  of  a  plane  wave  l>y  an  infinitely  staggered  guide 
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For  |fl|  >  tt/ 2,  $  satisfies  Helmholtz's  equation 


32$s  +  324>S  +  k2$s 
x  z 


where  k  ■  +  ik2,  k^  >>  k2  .  Introducing  the  Fourier  Transform 


(III. 3) 


4>S(x,a)  “  — 

/2ir 


,s,  N  iaz  , 
t>  (x, z)e  dz 


(III. 4) 


where  a  *  a  +  ir  and  <ts(x,a)  is  analytic  in  the  region  -k2  cos  JJ  >  x  >  -k,,, 
we  get 


32$S(x,a)  -  y2QS(x,a)  -  0 


(III. 5) 


2  2  1/2 

where  y  ■  [a  -  k  ]  is  defined  on  the  upper  sheet,  shown  in  Figure  3. 

The  solution  (III. 5)  satisfying  the  radiation  condition  can  formally 
be  written  as 


4>  (x,ci) 


x  >  0 


A(a)e"YX 

B(a)e  ^  +  C(a)eYX  -b  <  x  <  0 


(III. 6) 


From  the  continuity  of  4>  (x,a)  across  a  PEC  (Perfect  Electric  Conductor) 
at  x  ■  0,  one  obtains 


A(a)  =  B(a)  +  C(a) 


If  the  bottom  wall  is  a  PEC  as  well, 


(III. 7) 


-2yb 


B(a)  -  -C(a)e 


and  if  a  PMC  (Perfect  Magnetic  Conductor)  wall 


(III. 8) 


B(a)  -  C(a)e"2Yb 


(III. 8b) 
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If  J  (a)  Is  defined  as 


J_(a)  «  y{3x$S(0+,a)  -  3x$s(0“,a)}  -  -y(A  -  B  +  C)  , 


then  from  (III. 7),  (III. 8a)  and  (III. 8b)  the  relationship 


'»  '  ^  fjW)  A(0)>  (+. 


r 

I’PEC  wall' 

l^PMC  wall 

is  obtained.  If  the  boundary  condition  is  applied  at  x  *  0,  the 
Wiener-Hopf  Equation 

-4^(0, a)  +  4^(0, a)  *  J_(a),  -k2  cosfi  >  t  >  -k2 

is  constructed,  where 

G(o) 


and 


*i(0,o) 


1  -  e-2^b 

PEC  wall 

1  +  e'2yb  =  G(a) , 

PMC  wall 

-i(l  +  e_2ikbsinfl) 

2/2tT  (a  +  k  cos  0) 

t  < 

in  Section  A,  Chapter  II,  the  t 

(III. 9) 


(III. 10) 


(III. 11) 


(III. 12) 


(III. 13) 


process  is  the  factorization  of  G(a)/y.  The  factorization  of  y  is  obvious. 


and  from  Equation  (II. 4),  G(a) 

factors  into  G+(a)G_(o),  where 

(a)  For  G(a)  ■  1  -  e~2^b 

t 

(PEC  wall) 

1/2 

G+(a)  -  (2b  sinc(kb) \  4 

(a  + 

k)1/2  e"i1T/4  exp[~(l  -  C  +  In 

e»p(^ln(2-^)] 

00 

n 

— _  i 

(1  +  J-)  eiab/n\  r  >  -k2. 

—  +  — )  ] 
kb  2  n 

(III. 14a) 
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(b)  For  G(a)  -  1  +  e~2yb  =  G(a) ,  (PMC  wall) 


G+(a)  -  [2  cos  (kb)]^2  exp[^|^(l  -  C  +  +  -y-)] 


exp[“p-  ln(a  -"■■*)  H  (1  + 


a _ .  ictb/(n-l/2)ir 


-1  ^n-1/2 


)  e 


(III. 14b) 
,t  >  -k0 


(c) 


G_(a)  *  G+(-a) 


(III.14c) 


In  the  above  equations  C  is  the  Euler's  Constant  and  iYn>  iYn  are 
-2yb  ^2vb 

the  zeroes  of  (1  -  e  )  and  (1  +  e  ) ,  respectively,  with  yr  given  by 


r,mr,2  .2,1/2 
Yn  *  t(y-)  -  k  ] 


(III. 15) 


with  the  same  branch  cuts  as  in  Figure  3.  A  set  of  tables  and  other 
forms  for  G^(a)  are  presented  in  Appendix  B. 

Now  decomposing  Equation  (III. 11)  and  using  the  edge  condition,  the 
scattered  field  in  the  transform  domain  for  -b  <  x  <  0  is 


4>S  (x,a) 


(1  +  e~21kb3lnfl)  x(kcosfl.a) 


r^T 


2Y 


G+(-kcosft)G+(-ct) 


/.  t  -2ikbsin0. 
+  (1  ±  e  ) 


'2tt 


e-Y(xf2b) 

' 

1 

2y 

G  (-kcosft)G  (-a) 

l  +  + 

(;)■ 


PEC  wall 
PMC  wall 


(III. 16) 


where  X(k  cos  ft, a)  is  the  spectral  diffraction  coefficient  defined  in  (11.26). 
Comparing  Equations  (III. 16)  and  (11.27),  it  is  apparent  that  the  spectral 
representation  of  the  scattered  field  has  the  same  form  as  that  for  the 
isolated  half-plane  and  its  image,  but  the  half-plane  solutions  are  multiplied 
by  the  factor  in  the  curly  brackets,  which  accounts  for  the  multiple  inter¬ 
actions  between  the  top  plate  and  the  PEC  or  PMC  wall.  We  should  note  that  as 
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b  G+(a)  -*•  1  and,  since  Re(y)  >  0,  the  reflection  term  on  the  right-hand 

side  of  (III. 16)  goes  to  zero,  so  that  only  the  half-plane  solution  remains. 

To  find  the  field  in  the  guide,  the  inverse  Fourier  Transform  of 
(III. 16)  is  computed.  Since  the  guide  is  in  the  region  z  <  0,  the  integration 
path  may  be  closed  as  shown  in  Figure  12  and  the  residues  evaluated. 

In  computing  these  residues  one  finds  that,  first,  the  contribution 
around  the  branch  cut  is  zero.  Next,  the  residue  at  o  ■  -k  cos  ft  gives 
the  negative  incident  field  so  that  the  total  field  in  the  guide  is  fully 
expressible  in  terms  of  its  modes,  as  it  must  be.  Lastly,  the  poles  at 
a  *  iyn  give  the  modes  of  the  guide.  Note  that  a  ■  iy^  gives  the  TEM  mode. 

Summarizing  the  final  results,  if  the  field  of  the  guide  in  Figure  10 
is  written  as 

^(x.z)  -  l  C  sinfe^  x 

n-0  n  > 

where 

kg  -  [kg  -  [f1]2]172  ,  (III.  17bb) 


-ik  z 
e  g 


1.2, 


(III. 17a) 


then  the  modal  coefficient  Cn  is  given  by 


,  i(n+l)ir  -ikasinft-. 
(1  +  e  _ e _ 2. 

2ka  cos  $ 


D(fl,$n  -  ir)^ 


[G+(k  cos  $  ) 


G+(-k  cos  ft) 


(III. 18) 


with  n  even  coinciding  with  the  asymmetrical,  PEC  wall  problem  in  Figure  10a 
and  n  odd,  G+  *  with  the  PMC  wall  problem.  D(ft,4>n  -  ir)  is  the  diffraction 
coefficient  given  in  Equation  (11.48)  and  is  defined  by 


(III. 19) 


I 


«  cri- 


sin  *n  -  c 
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For  (niT/ka)  <  1,  that  is,  for  propagating  modes  we  can  interpret  (III. 18) 
to  be  the  amplitude  of  the  edge  diffracted  ray  in  the  modal  direction  of  the 
guide,  Figure  13,  where  2ka  cos  term  serves  as  a  normalization  factor. 
However,  instead  of  just  the  diffraction  coefficient,  the  "modified  diffraction 
coefficient"  [19] 

D(fl,9)  -  D(fi,9)f(p  ♦  nJfC-k1  •  n)  (III. 20a) 


where 


f(x) 


'G+(M> 

l/G+(x) 


x  <  0 
x  >  0 


(III. 20b) 


(see  Figure  14)  is  used  to  account  for  the  interactions  between  the 

plates.  Although  a  ray  picture  is  not  totally  proper,  we  say  that  the 

nature  of  this  interaction  is  the  reflection-diffraction  cycle  that  takes 

place  along  the  aperture,  as  shown  in  Figure  15.  This  interpretation 

is  not  strictly  valid,  because  after  the  first  reflection  off  the  top 

plate,  the  field  about  the  aperture  is  no  longer  a  ray  field  and  the 

reflection-diffraction  phenomenon  cannot  really  be  localized  to  a  point. 

To  summarize  we  note  that  (III.  18)  is  finite  even  when  D(ft,  -  tt) 

n 

becomes  unbounded  and  also  that  the  derivative  of  with  respect  to  a  is 
discontinuous  at  the  mode's  cutoff  point.  For  an  H-wave,  the  same 
considerations  are  used  but  the  diffraction  coefficient  is  that  for  the 
H  wave. 

B.  Modified  Diffraction  Coefficient  (MDC) :  Large  Finite  Stagger 

Having  examined  in  the  previous  section  the  unstaggered  case,  a  similar 
Wiener-Hopf  formulation  is  possible  for  the  finite  stagger  in  Figure  1. 

This  time,  however,  there  is  a  need  to  consider  an  additional  unknown,  the 
current  on  the  bottom  plate.  This  in  turn  leads  to  a  set  of  coupled 
Wiener-Hopf  equations  for  the  two  plate  currents  which,  in  general,  are  not 


»  T  -vr  Tj 


* 
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amenable  to  a  solution.  But  if  the  diagonal  distance  "d"  is  large  with 
respect  to  the  wavelength,  asymptotic  methods  may  be  employed  to  decouple 
the  Wiener-Hopf  equations,  and  an  approximate  solution  for  the  modal 
coefficients  can  be  obtained.  A  short  presentation  leading  to  the 
Wiener-Hopf  equations  is  now  given  and  the  asymptotic  solution  for  the 
modal  coefficients  presented.  The  details  of  this  section  can  be  found 
in  [19]. 

We  consider  again  the  incident  field  to  be  a  plane  E-wave  given 
by 

/(p,8)  -  eikpcos(n"6)  ,  (III. 21) 

with  the  same  restriction  -  |ft|  >  ir/2. 

The  scattered  field  is  defined  as  in  (III. 21),  and  with  the  Fourier 
transform  given  in  (III. 4)  one  again  arrives  at 

32$s(x,a)  -  y2$S(x,a)  ■  0  .  (III. 22) 

The  solution  to  (III. 22)  is 


A(ci)e  YX 

x  >  0 

<J>S(x,a)  *  ■ 

B(a)e-yx  +  C(a)eYX 

-a  <  x  <  0 

(III. 23) 

D(a)e+YX 

x  <  -a 

Continuity  of  the  E  field  gives 

A(a)  -  B(a) 

+  C(a) 

(III. 24) 

B(a)eya  +  C(a)e_Ya  -  D(a)e“Ya 

(III. 25) 
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The  top  plate  current  is  again  defined  as 

JT(ct)  -  J*(«)  -  Yt3x$S(0+’a)  "  3x*!(0“,ct)]  *  ~2  (A  -  B  +  C)  ,  (III. 26) 

where  the  superscripts  T  and  B  indicate  the  top  and  bottom  plates, 
respectively. 

Introducing  a  shifted  Fourier  transform 


<t>3(x,a)  =  e-iai  <J>3(x,a)  , 


(III. 27) 


we  observe  that 


*S(x,a) 


1 

/2ff 


,s,  .  ia(z-S,). 

(j>  (x,a)e  dz 


(III. 28) 


is  analytic  in  the  region  T  <  cos  ft.  Then  making  use  of  (III. 28),  the 
bottom  plate  current 


JB(a)  -  -|-[3x<t3(-a+,a)  -  3x$S (-a”, a)  ] 


(III. 29) 


can  be  written  as 


Aa)  -  e+la,?(a)  -  -  l[B?a  -  Ce-Ta  +  De-yal 


(III. 30) 


where  J  (a)  is  analytic  in  the  lower  half -plane  given  by  t  <  cos  £2. 

Two  additional  equations  are  obtained  by  applying  the  boundary 
condition  to  the  top  and  bottom  plates: 


-$^(0,a)  +  $3(0,a)  ”  A(a) 


(III. 31) 


and 


eiaS,[-$^(-a,a)  +  $3(-a,a)] 


D(a)e“Ya 


(III. 32) 
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s 

After  some  manipulation  and  use  of  the  edge  condition,  the  coupled 
Wiener-Hopf  equations 

<fra(0,a)  +  +  e"Ya  eia£J®(a)]  -  $*(0,cx)  -  0  (III. 33a) 

and  -k2  cos  ft  >  t  >  -k2  . 

eiaV^(-a,a)  +  i[e"yaJ^(a)  +  ei<UJ®(a)  ]  -  eiai^(-a,a)  -  0  (III. 33b) 

are  obtained.  It  is  obvious  now  that  for  an  arbitrary  stagger  i  the 
Wiener-Hopf  procedure  becomes  intractable.  However,  as  stated  before,  an 
asymptotic  evaluation  of  Equation  (III. 33)  can  be  made  and  the  field 
in  the  waveguide  computed  to  0(l/kd) .  The  details,  which  can  be  found  in 
Lee  [19],  are  quite  involved  and  only  the  final  result  is  now  given. 


If  the  field  in  the  guide  is  given  by 


f>C(x,z)  *  l  C  sin 


n-1 


mr 


-ik  z 
»  8 


n  *  1,2,3  . . . 


(III. 34) 
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r.2  fair)  2. 1/2 

[k0  '  (T“J  ] 


(III. 35) 


then,  with  H(A)  defined  as 
f 1  if  A  is  true 


H(A) 


0  if  A  is  false 


the  asymptotic  solution  for  for  the  propagating  modes  is 


C  -  I  cp}  +  0(l/kd) 

n  u  Ti 


n  j-i  n 


(III. 36) 


(III. 37) 


*  „  l  ■**-  sh*  r,  . 
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where 


.(1)  ,  [H(fl  <  -6)  -  e~2ikaslnfl  H(n  >  <5  -n)] 
'n  2ka  cos 


D(n,4>ti  -  it) 


G  (k  cos  <J>  ) 
+  n 


G  (-k  cos  Q) 


(2)  „  [H(*  <  6)eikaCOS(n+,5)][D(a,TT  -  *  )] 

n 


,  ,.n  ik£cos<{> 
(-1)  e  n 


2ka  cos  <j> 


(III. 38a) 

(III. 38b) 


.(3)  ,  felfcacoa(W-a)g(kd)D(ntir  _  &)] 


2ka  cos 


D(-6,<t>n  -  tt) 


G,  (k  cos  <p  ) 
“p  n 


G+(k  cos  6) 


(111.38c) 


and 


c(4)  -  H(<fr  <  <5)[H(fl  <  -6)-e“2ikasinil  H(«  >  6  -  tt)  ]  (g(kd)D(fl,  5  -  n)  ] 
n  n 


[  D  ( TT  -  6,  TT  -  (J»n)] 


f  1  1 

f(-l)Vik*cos*n] 

G+(-kcos  Q)G+(kcos 

5) 

t 

2kacos  $ 

n  J 

(III . 38d) 


The  angle  is  defined  as  in  Section  A  and  is  given  by 

sin  <t>  -  ~  •  (III. 39) 

n  ka 

For  propagating  modes  the  four  numbers  c^^  thru  c ^  can  be  given 
the  interpretation  of  different  diffracted  rays  "coupling"  into  the  guide. 
Consider,  for  example,  c^^  .  We  see  that  c^^  is  the  amplitude  of  the  top 
edge  diffracted  ray  in  the  modal  direction  -  the  same  interpretation  that 
was  given  in  Section  A.  The  functions  H(*)  in  the  equation  for  c^^  are 
present  so  that  a  ray  "couples"  into  the  guide  via  the  edge  diffraction 
only  if  that  ray  hits  the  edge.  The  situation  is  illustrated  in  Figure  16. 

We  note  the  presence  in  (III. 38)  of  the  "Modified  Diffraction  Coefficient," 
which  as  in  Section  A,  is  said  to  account  for  the  interaction  between  the 
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/ 


(2) 

plates  along  the  guide  aperture.  The  interpretation  of  cn  is  similar,  but 

ik&cosd) 

with  an  extra  phase  factor,  e  n,  present.  The  diffraction  coefficient 

is  now  unmodified,  as  the  bottom  plate's  edge  is  not  on  the  aperture.  The 

terms  c^  and  c^  are  due  to  the  twice-diffracted  diagonal  rays,  shown  in 
n  n 

Figure' 16.  Higher-order  interactions  along  the  diagonal  are  ignored,  as  this 
result  is  only  0(l/kd). 

In  Sections  A  and  B,  we  have  presented  the  modal  coefficients  and 
have  associated  them  with  the  edge  diffraction  rays.  But  the  mechanism 
by  which  the  diffracted  rays  couple  into  the  guide  and  the  reason  for 
the  absence  of  the  incident  field  in  the  final  solution  are  not  clear. 

This  and  other  points  are  discussed  in  Section  C  and  D. 

C.  YKF  Method 

Though  the  exact  solution  of  the  unstaggered  case  was  given  in 
Section  A,  a  ray  optics  solution  to  the  problem  is  nonetheless  highly 
desirable.  Not  only  are  we  provided  with  another  check  of  the  different 
asymptotic  theories,  but  the  use  of  ray  optics  and  the  local  nature  of 
ray  methods  allow  for  a  better  physical  understanding  of  how  the  energy 
is  actually  coupled  into  the  guide.  The  first  published  work  along  this 
line  was  by  Yee  and  Felsen  [ 30 ] >  and  in  the  ray  analysis  presented 
below,  we  often  refer  to  this  paper. 

C.l.  Formulation  via  ray  optics 

We  start  by  considering  the  field  in  the  presence  of  the  unstaggered 
guide.  For  exactness,  the  incident  field  is  an  E-polarized  line  source 
field,  normalized  so  that  its  phase  center  is  at  (0,0).  That  is, 


4>i(P„)  -  6 


ik(pg-R) 


(III. 40) 
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where  p  is  the  distance  from  the  source,  as  in  Figure  17.  The  advantage 
s 

of  using  a  line  source  in  our  formulation  becomes  evident  later. 

If  R  is  large  enough,  we  have  in  essence  a  plane  wave.  Using 
Equation  (11.45),  the  total  field  outside  the  guide  may  be  written  as 

/  =  [F(5^)  -  F(£*)  +  F^)  -  F(£*)  “  !!♦* 


-  [F(eJ)  -  FaJmJ  -  [F(C  J)  -  Fa^)  ]<f>2 

+  ($J  +  $jnt)  +  (**  +  $2nt)  •  (III. 41) 


In  the  above  equation,  the  subscripts  1  and  2  refer  to  the  plates  as 
given  in  Figure  17,  £*’!l  are  the  detour  parameters,  7  are  the 
diffracted  fields  that  would  exist  if  the  plates  were  infinitely  far 
apart,  and  '<>*^2  are  tile  Cerms  due  co  che  interaction  between  plates  1 


We  call  the  GO  and  diffracted  fields,  without  any  interaction, 
the  "primary  field."  The  -$*  term  in  the  first  bracket  in  (III. 41)  is 
present  so  that  the  incident  field  is  not  counted  twice.  We  should  emphasize 
that  Equation  (III. 41)  is  only  for  the  region  outside  the  guide.  To  find 
the  field  inside,  the  guiding  action  of  the  waveguide,  that  is,  the  reflec¬ 
tions  of  the  field,  as  schematically  indicated  in  Figure  18,  must  be  consid¬ 
ered.  Once  we  have  done  so,  and  the  far  field  inside  the  guide  is  known, 
the  model  coefficients  are  computed  by  simply  Fourier  decomposing  the  field 
over  the  aperture  A.  That  is 

-  ra  ikzcos<$ 

C  »  —  <f>  (x,z)  sin(ksin  $  x)e  n  dx  (III. 42) 

n  3.  j  n 


We  have  illustrated  this  procedure  by  first  considering  the  diffracted 
fields  and 


-  - 


V- 


gap** 


9 


43 

C.2.  Primary  field;  Diffracted  field  contribution 

Consider  the  line  source  at  x  *  0,  z  ■  0,  producing  in  free  space 
the  diffracted  field 

^  -  ^(O,  O)D(fl,0)g(k p);  (III. 43) 

g  and  D  are  defined  in  Equations  (11.47)  and  (11.48).  If  one  ignores  any 
diffraction  by  the  lower  plate  and  just  traces  the  reflected  rays  in  the 
guide,  the  far  field  in  the  guide  excited  by  the  edge  line  source  is  the 
same  as  that  produced  by  the  infinite  set  of  images  located  at  z  *  0, 

x  *  2a j  j  -  0,  ±1,  ±2  .  (III. 44) 

That  is,  for  z  <  0,  -a  <  x  <  0,  the  edge  line  source  gives  rise  to 

a  far  field  given  by  the  field  of  an  infinite  linear  array  centered  at 

the  edge  and  spaced  2a  apart,  (see  Figure  19),  with  the  array's  field  given  by 

03 

♦2  -  ^(0,0)  l  D(fl,0  )g(kp  )  (III. 45) 

j”-®  J  J 

where 

re  o  ^2  ,  2,1/2 

Pj  “  [(x  -  2aj)  +  z  ] 

and 

*  tan-1[(x  -  2aj)/z]  ,  j  9^  |  <  tt  . 

Note  that  since  D(ft,  -  9^)  «  -D(fi,  0^)  the  array  is  antisymmetrical.  With 
the  modes  as  given  in  (III. 17),  the  contribution  to  the  modal  coefficient 
from  <}>2  is 
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o  *  i  ra  ^  d 

r  -  —  4>“(xfz)  sin(k  sin  6  x) 

n  a  j  l  n 


ikzcos<$> 


(III. 46) 


with  z  »  0.  If  the  change  of  variables  t  ■  x  -  2a j  is  introduced, the 
integration-summation  in  (III. 46)  becomes  an  infinite  integral  and 


..  ,  .  ikzcos<(i 

ikp(t)  .  ,,  .  ,  .  n 

e  sin(k  sin  $  t)e 


dt.  (III. 47) 


If  the  integral  is  evaluated  asymptotically,  we  obtain 

,  «i(0,0)D(n,  A  -  it) 

CP  -  - S - 

n  2ka  cos  <t> 


(III. 48) 


Likewise,  for  the  line  source  at  x  *  -a,  we  get 
o'  i(n+l)ir.i,n  .  D(n>  ^n  ~ 

cn  e  *  (0*  "a)  “2kT^7T  • 

n 


(III. 49) 


Adding  (III. 49)  to  (III. 48),  we  see  that  the  part  of  the  modal  coefficient  due 
to  diffraction  is  the  same  as  the  total,  exact  modal  coefficient  in  (III. 18) 
except  for  the  factor  in  (III. 18)  which  accounts  for  the  interaction  between 
the  plates.  Note  that  the  results  derived  here  are  only  valid  for  the 
propagating  modes,  whereas  (III. 18)  applies  to  the  evanescent  modes  as  well. 

C. 3.  Primary  field:  GO  field  contribution 

With  the  interpretation  just  given  for  Equation  (III. 18),  it  appears 
that  the  geometrical  optics  field  plays  no  role  in  the  final  solution. 

To  see  why,  we  set  up  the  problem  as  before.  First,  looking  at  (III. 41)  and 
Figure  17,  it  is  evident  that  for  the  region  in  the  waveguide  it  is  sufficient 
to  consider  the  geometrical  optics  field  given  by 


^  ^  ,  w  ry- 

..  '*  ' 


Vr  -  CT--  -- 


r 


$G  -  m*)*1  -  FUr2)tT2 


where  ,  for  convenience 


F(x)  =  F(x)  -  F(x)  . 
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(III. 50) 


(III. 51) 


That  is,  the  field  is  that  of  two  line  sources,  one  with  a  pattern  given 

i  r 

by  F(C^)  and  the  other  with  a  pattern  where  from  Figure  20 


4  -  eX’l^uJ  +  bJ-cJ) 
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(III. 52a) 


(III. 52b) 


(III.52c) 
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*5 

<  i 

i«2i 

(III.52d) 

-1 
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> 

laj 

0 

i  2 1 

and  6q 

are 

of 

opposite  sign. 

Note  that  6^ 

We  first  examine  the  source  F(5^)  and  look  at  its  array  of  images  for 
the  field  in  the  waveguide.  This  is  shown  in  Figure  21.  The  source  is 
at  z  »  -R  cos(01),  x  *  R  sin(|ftjJ);  the  Images  are  located  at 
Zj  *  -R  cos(O^),  Xj  *  -R  sin(|f2^|)  +  2aj ,  j  ■  ±1,  ±2,  ...  .  Observe 
that  for  x  <  0  the  pattern  is  -F(5^(-9j))  because  the  bottom  images 
account  for  an  odd  number  of  reflections.  This  array  is  allowed  to 
radiate  in  free  space  with  the  detour  parameters  given  by  (III. 52a) 
and  (III. 52b),  the  subscript  on  the  variables  going  from  0  to  j . 

A  similar  array  for  the  source  ~F(,Z^),  located  at 
z  ■  R  cos(01),  x  *  -R  sin(|n^|)  -  2a,  is  shown  in  Figure  22.  Again, 
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Figure  22.  Array  of  inage  sources  for  the  GO  field  F(57) 
of  a  source  at  x  ■  R  s in ( 1 2  j ) . 
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note  the  change  of  sign  for  x  >  0  and  also  that  now  there  is  no 

image  at  x  *  R  sin  (|o^|).  As  before,  the  detour  parameters  are  calculated 

from  (III. 52c)  and  (III.52d). 

Now  we  let  R  be  large  but  still  finite.  It  immediately  follows  that 
^1  ”  ^2*  furthermore,  if  we  let  the  observation  point  recede  deep  into 
the  waveguide,  that  is,  |z|  both  and  become  large,  but  as  can 

be  seen  from  (III. 52),  are  of  opposite  sign.  The  addition  then  of 
the  fields  in  Figure  21  and  22  leaves  the  field  from  the  source  at 
z  *  R  cos(JJ^),  x  *  R  sin(jftjJ),  radiating  in  free  space.  But  since 
is  now  large  and  the  observation  point  is  deep  inside  the  waveguide,  in  the 
shadow  region,  we  conclude  that  the  GO  field  deep  in  the  guide  is  zero. 

The  fact  that  neither  Equation  (III. 18)  or  (III. 48)  shows  the  GO  field 
should  not,  however,  be  taken  to  mean  that  only  the  diffracted  field  is 
needed  when  considering  the  primary  field  contribution.  We  can,  instead 
of  looking  into  the  far  field  of  the  guide,  find  cP  by  integrating  the 
primary  field  over  the  guide  aperture  at  z  »  0.  That  is. 


2 

a 


((F(^)  +  F(e*)  -  l)*1  -  F(eJ)^  -  F^*)*^ 


sin 


z-0 


mr 


dx 

(III. 53) 


where  we  used  (II. 34)  to  express  the  primary  field.  If,  however,  instead 
of  the  full  expression  in  (III. 53)  we  were  to  use  the  diffracted  fields 
only  and  calculate  cp  by 


p'  2 

*  r  tm  — 


-a 


sin 

[nr 

—  X 

z«0 

|a 

dx 


(III. 54) 


cP  would  be  greatly  underestimated.  This  can  be  seen  in  Figure  23,  where 


for  a  typical  case,  we  compare  the  magnitudes  of  cp  as  calculated  from  (III. 48) 
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and  (III. 49),  (III. 53)  and  (III. 54).  We  also  observe  that  the  asymptotic 
solution  breaks  down  near  cut-off. 

The  answer,  then,  to  the  question  of  whether  one  should  include 
the  GO  field  is  dependent  on  where  the  observation  point  is.  If  one  is 
examining  the  aperture  field,  the  GO  field  must  be  included.  If  one  is 
studying  the  far  field,  the  GO  field  drops  out,  hence,  the  appearance 
of  only  the  diffraction  coefficient.  But  this  does  not  mean  that  only 
the  knowledge  of  the  diffracted  field  is  sufficient  to  find  c^.  It  is 
sufficient  only  if  the  integration  is  done  in  the  far  field. 

C.4.  Primary  field:  Aperture  field  method 

Equation  (III. 53)  suggests  another  way  of  computing  the  primary 
contribution  -  find  and  integrate  the  primary  field  over  the  aperture. 

This  approach  has  several  important  advantages.  First,  by  not  looking 
down  the  guide,  there's  no  need  to  track  the  reflections  which,  as  just 
seen,  can  lead  to  complicated  forms.  And,  although  analytic  expressions 
such  as  (III. 48)  and  (III. 49)  cannot  be  found,  the  aperture  field  integral 
can  be  easily  and  very  accurately  numerically  integrated,  as  its  limits 
are  from  -a  to  0.  Furthermore,  though  Figure  23  seems  to  indicate  little 
reason  for  choosing  the  aperture  field  method  over  the  asymptotic  solution, 
the  aperture  field  method  is  superior  in  that  Equation  (III. 53)  contains  no 
singularities  and  hence  the  procedure  is  easily  extended  to  the  staggered 
case.  Additional  details  on  the  aperture  field  method  are  discussed  in 
Chapter  IV. 
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C.5.  Multiple  scattering 

The  solutions  in  Equations  (III. 48),  (III. 49)  and  (III. 53)  are  based 
on  the  assumption  that  each  plate  in  the  guide  is  excited  by  the  incident 
field  alone.  In  fact,  each  plate  is  also  excited  by  the  fields  scattered 
off  the  other  plate  so  that  the  effects  of  multiple  scattering  must  be 
considered.  In  the  exact  solution,  this  effect  has  been  accounted  for 
in  toto  by  the  use  of  the  Modified  Diffraction  Coefficient  (MDC) . 

The  ray  optics  solution,  on  the  other  hand,  given  by  Yee  and  Felsen  [30], 
attempts  to  keep  track  of  each  of  the  interactions.  This  is  now 
briefly  outlined  and  the  limitations  of  the  technique  discussed. 

Consider  for  example  the  field  scattered  off  the  top  plate,  as 
in  Figure  24.  We  denote  the  scattered  field  off  plate  1  after  (m  +  1) 
scatterings  by  <t>®  ,  likewise  for  plate  2  by  d>~  .  To  find  the  scattering 

■L  y  Ul  Z  y  Ul 

g 

of  q  by  plate  2,  the  YKF  solution  proceeds  as  follows.  First,  assume 

g 

that  the  initial  scattered  field  ^  q  near  the  lower  edge  is  well- 

approximated  by  the  diffracted  field.  Next,  approximate  this  first 

diffracted  field  by  an  isotropic  line  source  field,  its  strength  given  by 

D(fl,  ir/2).  As  this  field  scatters  off  plate  2,  the  resultant  scattered 

field  near  the  aperture,  denoted  by  ,  is 

Z9  X 

$2 *  DO*.  w/2){g(ka)D(*/2,92)g(kp2)  -  [F(^)  -  F(^)  ]g(kPr)  }.  (III. 55) 

As  we  again  consider  the  scattering  of  off  plate  1  to  obtain  - ,  we 

Z9  1  Ly  X 

observe  that  the  top  edge  and,  indeed,  the  entire  aperture  lies  on  the 

reflected  shadow  boundary  of  $9  ..  However,  note  that  even  on  the  shadow 

z ,  X 

s  1 

boundary  $2  ^  is  finite  and  is  equal  to  -^(ft,  ir/2)g(kpr). 

To  find  the  scattered  field  again  the  approximation  is  made  that 

<t>2  ^  is  an  isotropic  line  source  field  of  strength  -|d(^,  tt/2)  and  the  entire 

procedure  repeated  to  find  the  next  higher-order  scattering.  When  this 
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is  all  done,  the  modal  coefficient  due  to  the  interactions  is 

*  ■  '  -  r 

l 


c^nt  -  D(Q,  tt/2)  +  e1(n+1)%i(0,  -a)]  x 

n  m-1 


ima 


>m-l 


l*G  (2)" 

(III. 56) 


That  is,  the  multiply  diffracted  fields  are  reduced  by  one-half  upon  each 

diffraction  and  the  coupling  of  these  fields  computed  as  in  Section  C.2. 

An  interesting  comparison  of  (III. 56)  with  the  asymptotic  expansion  of 

the  exact  solution,  made  by  Bowman  [9  ],  shows  that  (III. 56)  would  tend 

to  underestimate  c*nt. 

n 

The  YKF  method  for  the  interactions  then  makes  use  of  three  basic 
assumptions: 

(a)  On  opposite  plates,  the  primary  scattered  field  is  well- 

approximated  by  the  diffracted  field.  This  will  be  true  if 

a  is  large  and  <j>  is  not  close  to  w/2. 
n 

(b)  The  diffracted  fields  are  approximated  by  an  isotropic  line 
source  field. 

(c)  A  greater  approximation  than  the  isotropy  of  the  diffracted 

fields  is  that  ray  fields  are  used  at  all.  Equation  (III. 55) 

clearly  shows  that  after  the  first  interactions  the  edges  lie 

in  the  transition  region  of  the  interacting  fields,  where  a 

ray  picture  of  any  kind  is  inappropriate.  Furthermore, 

near  cutoff,  that  is,  <j>  -*■  r/2,  the  ray  picture  becomes 

n 

even  more  confusing.  In  Chapter  IV,  the  interactions 
are  partially  taken  into  account  in  a  more  exact  manner. 

D.  Summary  and  Conclusions 

D*l*  Exact  and  asymptotic  Wiener-Hoof  solution  (MDC) 

In  Section  A  an  exact  solution  for  the  modal  coefficients  of  an 
unstaggered  parallel-plate  guide  was  found.  The  expression,  given  in 


i  v- 


£3*: 
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Equation  (III. 18),  shows  that  the  modal  coefficients  are  proportional  to 
the  amplitude  of  the  diffraction  coefficient  in  the  modal  direction, 
multiplied  by  a  factor  which  accounts  for  the  interaction  between  the 
plates.  This  then  leads  to  the  introduction  of  the  Modified  Diffraction 
Coefficient  (MDC)  in  Equation  (III. 20).  Note  that,  generally  speaking, 
the  MDC  will  not  give  a  ray  field  since  it  is  a  function  of  k,  indicating 
that  the  interaction  is  of  a  complicated,  nonray-like  nature. 

In  Section  B,  an  asymptotic  solution  valid  for  large  was 
introduced  in  Equation  (III. 38).  We  note  again  the  use  of  the  MDC, 
but  now  also  the  presence  of  singularities  in  Eouations  (III. 38c)  and 
(III.38d)  as  6  -*■  4>n.  These  singularities  arise  from  the  assumptions  of 
the  asymptotic  evaluations  and  may  in  principle  be  removed  [1],  though 
this  would  lead  to  exceedingly  complex  forms. 

Taking  into  account  the  inherent  limitations  of  the  asymptotic 
solution,  when  can  Equation  (III. 38)  be  expected  to  give  good  results? 
"When  "l"  is  large  enough  so  that  the  lower  edge  is  outside  the  transition 
region  centered  about  the  aperture,  as  shown  in  Figure  25."  This 
transition  region  is  the  region  in  which  the  interactions  between  the 
plates  take  place,  as  outside  of  this  region  we  have  ray  fields  which 
only  produce  specular  reflections.  When  the  lower  edge  falls  out  of 
this  region,  the  use  of  the  MDC  to  account  for  the  interactions  is 
justified.  Simple  calculations  show  "£"  would  have  to  be  greater  than 
a  wavelength,  for  example,  when  a  *  .55  X,  l  >  1.1  X.  When  a  =  1.2  X, 
the  minimum  "i,"  increases  to  l  >  1.4  X. 
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D. 2.  YKF  and  aperture  field  method 

In  Section  C,  a  ray  optics  solution  for  the  unstaggered  guide  was 
presented.  It  was  shown  that  the  primary  field's  contribution  to  the 
modal  coefficient  was  through  the  diffracted  fields  only,  the  GO  field 
being  absent  because  a  far-field  evaluation  was  involved.  It  was  also 
seen  that  because  the  guide  was  unstaggered  the  singularities  in  the 
diffraction  coefficient  cancelled.  An  equivalent  formulation  for  the 
primary  field's  contribution,  an  aperture  field  method,  was  given, 
and  shown  to  be  preferable  in  many  ways.  The  YKF  method  for  the 
interactions  was  discussed  and  its  shortcomings  pointed  out,  Tamely, 
the  presence  of  a  singularity  near  cutoff,  and  the  use  of  ray  fields  in 
the  transition  region. 

D.3.  Extension  to  the  staggered  case 

For  "i."  >  0  the  results  in  Section  B  have  an  obvious  and  immediate 
application.  However,  as  pointed  out  the  MDC  has  several  shortfalls. 
First,  for  the  MDC  to  be  successfully  applied,  the  stagger  must  be  large  - 
greater  than  one  wavelength.  Second,  the  terms  which  give  the  first 
interaction  along  the  diagonal,  Equations  (III, 38c)  and  (III.38d),  are 
singular  as  6  -►  it  -  fi  and  are  discontinuous  due  to  the  presence  of  the 
function  H(  ).  The  YKF  method  has  similar  problems;  a  singularity  when 
5  -*■  ir  -  fi,  the  need  for  "J,"  to  be  large  so  that  the  interaction  rays  are 
indeed  rays  on  the  opposite  plates,  and,  in  addition,  the  primary  field 
term  is  singular  when  $  -*■  it/ 2 .  A  comparison  between  the  MDC,  the  YKF, 

and  the  more  rigorous  Aperture  Field  Method  for  a  staggered  guide  is  given 
in  Section  D,  Chapter  IV. 
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IV.  STAGGERED  GUIDE 


We  examine  now  the  staggered  wave  guide,  l  +  0,  in  Figure  1.  As  in 
Chapter  III  we  only  study  the  E-wave  problem  in  detail  and  show  how  to  apply 
the  method  to  an  H-wave  case.  In  Section  A  we  find  the  plane  wave  primary 
field  contribution  to  the  modal  coefficient  and  make  use  of  the  Aperture 
Field  Method  introduced  in  Chapter  III.  In  Section  B  interactions  between 
the  plates  are  considered  and  the  Spectral  Theory  of  Diffraction  is  utilized 
to  exactly  solve  for  the  first  interaction.  Data  are  presented  in  Sections  C 
and  D  and  the  cylindrical  wave  problem  is  examined  in  Section  E.  Finally, 
we  note  that,  although  in  this  chapter  we  are  solely  interested  in  the 
propagating  modes,  the  same  methods  may  also  be  used  for  the  evanescent  modes 
Rather  than  dealing  with  C^,  the  modal  coefficient  of  the  n^  propa¬ 
gating  mode,  we  use  Pr,  the  power  transmission  coefficient,  defined  as  the 
ratio  of  the  mode  power  to  the  source  power  incident  on  the  aperture, 
projected  in  the  direction  of  the  waveguide,  i.e., 

C  | 2  cos 

P  „ — S - 5  ,  (IV. 1) 

n  2 j  cos (ft) j 

as  illustrated  in  Figure  26.  The  same  definition  is  used  when  ft  >  0,  though 
in  this  case,  has  less  physical  meaning.  Observe  that  because  of  the 
way  Pr  Is  defined  it  is  possible  for  it  to  be  greater  than  one.  We  define 
the  phase  shift  of  the  mode  as  the  phase  of  Cn  minus  the  incident  field's 
phase  at  the  center  of  the  aperture,  (-a/2,0).  That  is, 


\  a 

Phase  Shift  *  arg  -  arg  $(-■£,  0). 


(IV. 2) 


A.  Primary  Contribution  (Plane  Wave):  Aperture  Field 

The  primary  contribution  to  the  modal  coefficient  is  found  by  the 


Fourier  decomposition  of  the  primary  field  over  the  guide  aperture,  as 
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in  Equation  (III. 53).  But  now,  for  any  l  _>  0,  we  have  for  the  primary  field 
contribution 


cj  «  |  [  £(f<-4>  +  “  f(5F)$£  -  sin 

J -a  1 z*0 


1  sln 

nir 

—  X 
a 

*  z*0 

dx 

(IV.3) 


where 


^’r(x,0)  -  e^,r 


^2k)x|"  sin^j  +  -j 


C^’r(x,0)  «  E 


i,r 


Ak[(a  +  x)2  +  I1]111  cos  -  +  x  tan"1 


2-2 


a  +  x 


4>^(x,0)  -  exp[-kx  sin  Q] 


^(xt0)  *  exp[-ik(x  +  2a)  sin  0] 

When  l  -  0,  (IV.3)  reduces  to  (III. 53). 

In  Figures  27  and  28,  we  show  for  4*0  some  typical  comparisons  of 
Pr  and  the  phase  shift,  as  calculated  from  (IV.3),  and  the  same  parameters 
obtained  from  the  exact  solution.  We  see  that  (IV.3)  gives  accurate  results, 
the  error  in  Pfl  is  kept  under  10%  or  +0.8  dB,  and  the  phase  shift  less 
than  10°.  The  greatest  errors  occur  at  the  onset  of  a  mode  when  P 

n 

is  most  sensitive  to  small  changes  in  the  aperture  field.  As  a/X  increases, 
the  agreement  with  the  exact  solution  becomes  excellent. 

B.  First  Interaction  (Plane  Wave) 

The  fields  on  the  aperture  due  to  the  first  interaction  can  be 
solved  for  exactly  using  STD.  The  scattered  field  of  the  first  inter¬ 
actions  can  be  written  as 

»illt  "  *1,1  +  *2,1  <1V*4> 
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Figure  28.  Power  transmission  coefficient  and  the  phase  shi 
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s 

where,  following  the  notation  in  Section  C,  Chapter  III,  ft^  ^  is  the  field 

s  s 

scattered  off  plate  1  when  the  incident  field  is  q,  and  ^  is  the  field 

scattered  off  plate  2  when  the  incident  field  is  The  scattered  fields 

1,  u 

s  s 

q  and  $2  q  are  the  primary  scattered  fields  off  plates  1  and  2,  respectively, 
as  indicated  in  Figure  29.  From  Equation  (11.28)  the  primary  scattered 
fields  may  be  written  as 


L  S 

*1(2) 


ft)e 


ikpl(2)COS^“  i  ®i  I  ^ 


1(2) 


dft 


(IV. 5) 


where  (p^,6^),  (p2»02)»  an<^  l®l  <  ^  are  c^e  cylindrical  coordinates  attached 
to  the  edges  of  plates  1  and  2,  respectively,  and  r  is  the  path  in  Figure  5. 
Note  that  Equation  (11.28)  applies  when  ft  >  0.  For  ft  <  0  we  use  the  fact 
that  <ps  (—  |  ft  | ,  0)  equals  ij>s(|ft|  ,-0) ,  which  for  an  E-wave  is  the  same  as  <f>s(|ft|,0). 
Next,  using  superposition,  as  in  (11.39),  we  obtain  the  expressions  for 


>i,i  and  *5fi 


bs  _  ft, 

Pl,l  4iri 


X(|ft|,4-)e'ikd  C0S  (^5){F(-^)eikpl  C0S 


ikp  cos  (ft+0.) 
+  F(?*)e  1  i 


dft  , 


x  >  -a 


(IV  .6) 


where 


*  +>/2kp1  sin 


"ft  +  e. 


and 


*2,1  *4^  {  x(N.ft)e 


■  ikp„cos(ft+0„) 
ikd  cos  <*-«)|F(_;i)e  2  2 


ikp,  cos  (ft-0,) 
+  F(5')  e  l 


dft  ,  x  >  -a 


where 


$2,T  “  +^2kP2  sin 


ft  +  02 


(IV. 7) 
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Figures  29a  and  29b.  First  interaction  between  che  plates 


for  the  integral  in  (IV. 7).  The  saddle  points  associated  with  the  trans¬ 
formations  (IV.  9)  and  (IV. 10)  are  ir  -  6  and  6,  respectively ,  and  the 
residues  to  be  considered  are  from  the  poles  at  ^  »  +f2,  respectively. 

Using  then  the  steepest  descent  path  transformations  and  the  program  for 
complex  Fresnel  functions  in  Appendix  A,  c^nt  can  be  computed  with  any  of 
the  standard  quadrature  methods. 

The  effects  of  the  first  interaction  are  ex&si-'ned  for  the  test  case 

1*0.  The  same  examples  as  in  Figures  27  and  28  are  now  presented  in 

Figures  30  and  31  but  with  the  first  interaction  included.  We  see  that  the 

maximum  error  in  P  has  been  reduced  by  a  factor  of  two  to  under  0.3  dB, 
n 

and  under  5°  for  the  phase  shift.  Also  the  convergence  to  the  exact 
solution  as  a/X  increases  is  more  rapid,  with  the  agreement  already 
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Figure  30.  Power  transmission  coefficient  and  the  phase  shi 
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excellent  for  a/X  -  0.6.  Higher-order  interactions  are  hereafter  ignored  as 
an  exact  formulation  of  them,  possibly  via  an  integral  equation,  would  little 
improve  accuracy  and  any  approximate  technique  would  have  an  error  at  least 
as  large  as  the  error  present  now.  We  emphasize  that  the  formulations  here 
for  the  primary  and  interaction  terms  are  equally  valid  for  any  "fc"  and  the 
same  order  of  accuracy  is  expected  as  that  for  l  =  0. 

C.  Numerical  Results 

We  study  the  staggered  guide  making  use  of  (IV. 3)  and  (1V.8).  Figures  32 
thru  38  present  for  different  angles  of  incidence  ft  the  power  transmission 
coefficient  and  the  phase  shift  for  a  =  0.7  X,  0.8  X,  and  0.9  X  as  a  function 
of  the  stagger  distance  These  cases,  0.5  X  <  a  <  1  X,  are  of  particular 

interest  as  guides  in  the  dominant  mode  are  often  used  as  the  elements  in 
lens  antennas  [28]  and  in  other  antenna  systems  [14],  Furthermore,  it  is  for 
0.5  X  <  a  <  1  X  that  a  more  exact  analysis  is  needed;  for  plate  spacing  this 
small  asymptotic,  high  frequency  methods  cannot  be  expected  to  give  very 
good  results. 

Figure  32  presents  data  for  ft  »  180°,  that  is,  the  wave  vector  is 
parallel  to  the  guide  aperture.  As  2,/X  increases,  we  see  gradually 
decreasing,  falling  at  1  *  1  X  by  about  3  dB.  This  is  because  the  scattered 
field  on  the  aperture  is,  essentially,  the  diffracted  field  from  the  edges 
of  the  two  plates.  As  one  of  the  edges  moves  away  from  the  aperture,  the 
scattered  field  on  the  aperture  becomes  reduced.  However,  as  |ft|  decreases 
and  the  wave  vector  tends  to  parallelism  with  the  aperture,  several  different 
effects  can  be  discerned.  First,  for  l  *  0,  decreases  as  |ft|  decreases. 
This  is  because  as  |ft|  decreases  from  180°,  the  wave  number  k  sin  ft  increases 
in  magnitude,  which  results  in  a  more  rapid  phase  variation  of  the  incident 
field  on  the  aperture.  In  turn,  this  implies  less  power  coupling  into  the 


fundamental  mode. 


1.0 


‘  '  * 


Second,  for  ft  <  0,  that  is,  the  wave  incident  from  above,  as  |ft|  decreases 
and  1/ A  increases,  we  see  more  power  in  the  dominant  mode.  Although  the 
diffracted  field  on  the  aperture  from  the  bottom  plate  decreases,  more  energy 
is  now  channelled  into  the  guide  by  the  stagger  via  reflection.  On  the  other 
hand,  for  ft  >  0,  less  field  appears  on  the  aperture  as  the  stagger  now  acts 
to  block  out  the  incident  field.  For  |ft|  close  to  90°,  however,  this  pattern 
breaks  down.  The  behavior  of  P^  and  the  phase  shift  becomes  more  erratic  as 
the  entire  aperture  now  lies  in  the  transition  region  of  the  incident  and 
reflected  fields.  Note  that  as  |ft|  -*•  90°  the  projected  aperture  goes  to 
zero  and  hence,  from  (IV. 1),  P^  goes  to  infinity. 

Perhaps  the  most  interesting  result  is  the  phase  shift  behavior.  Except 
for  |ft|  close  to  90°,  we  see  that  the  phase  shift  is  very  weakly  dependent  on 
and  is  less  than  20°.  One  could  essentially  compute  the  phase  shift 
for  l  m  0,  for  which  an  exact  solution  is  available,  and  use  the  result  for 
l  i  0.  Tables  that  can  be  used  to  compute  the  phase  shift  for  l  »  0  are 
provided  in  Appendix  B.  It  is  because  of  this  very  weak  behavior  of  the 
phase  shift  with  "1"  that  simple  phase  front  designs  for  waveguide  lens 
antennas  are  successful  [28], 

D.  Comparison  of  the  Aperture  Field.  MDC  and  YKF  Method 

We  now  compare  the  aperture  field  method,  the  MDC  method  as  given  by 
Equation  (III. 38),  and  the  YKF  method.  The  YKF  formulation  for  the 
staggered  guide  is  the  same  as  the  MDC  formulation  in  (III. 38),  but  instead 
of  the  G+  functions,  Equation  (III. 56)  is  used  for  the  interactions  along 
the  aperture.  The  results  are  shown  in  Figures  39  and  40.  We  observe 
that  for  small  and  moderate  the  MDC  and  YKF  methods  fail.  This  is 
because  they  include  the  interaction  rays  along  the  diagonal  which,  as  1 
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Figure  ^0.  Magnitude  and  phase  shift  of  C, ,  ■  - 
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approaches  zero  must  be  somehow  removed,  as  they  are  not  present  in  the 
canonical  solution  for  A  ■  0.  We  observe  in  the  figures  the  singularity 
at  6  *  it  -  |fl|,  something  not  present  in  the  aperture  field  method.  As 
A  increases  all  the  methods  approach  the  proper  solution.  It  would  appear 
from  Figures  39  and  40  that  if  one  is  away  from  the  singularity  at 
5  »  it  -  |fl|,  the  MDC  and  YKF  methods  give  good  results  for  A  >  2X. 

E.  Incident  Wave-Source  Field 

An  advantage  of  the  aperture  field  method  is  that  it  allows  one  to 

consider  incident  fields  other  than  a  plane  wave.  As  an  example  we 

consider  the  problem  of  an  incident  line-source  field,  the  geometry  of 

which  is  given  in  Figure  17.  The  guide  is  unstaggered  and  the  source 

coordinates  are  (R  cos  £2^,  R  sin|ftj),  where  R  is  now  a  finite  number  of 

wavelengths.  The  problem  of  finding  the  primary  field  contribution  to 

the  modal  coefficient  can  be  attacked  in  one  of  two  ways.  One  could 

proceed  in  the  manner  of  the  YKF  solution,  as  in  Section  C,  Chapter  III, 

but  in  this  case  an  additional  singularity  is  introduced  in  Equations 

(III. 48)  and  (III. 49)  when  £2_  or  £2-  equals  <J>  -  ir.  We  recall  that  for  an 

1  2  n 

incident  plane-wave  the  singularities  in  (III. 48)  and  (III. 49)  cancelled 

out  whereas  for  a  line-source  field  they  do  not. 

On  the  other  hand,  the  aperture  field  method  is  simple.  Though  an 

exact  solution  to  the  problem  of  scattering  of  a  line-source  field  is 

available  [8],  we  prefer  to  use  the  asymptotic  solution  given  in  [7].  We 

begin  by  expanding  the  line-source  field  in  an  asymptotic  series.  Since  we 

-3/2 

utilize  the  primary  fields  only  our  results  will  be  accurate  up  to  0(k  ), 

-1/2 

and  hence  only  the  k  term  in  the  incident  field  expansion  is  retained. 

A  higher  order  of  accuracy  can  be  obtained  by  first  taking  the  nlane-wave 

spectrum  of  the  incident  field  and  for  each  plane-wave  component  computing 

*  ■* 
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Che  interactions  between  the  plates  as  in  Equations  (XV. 6)  and  (IV. 7).  This, 
however,  would  prove  to  be  too  laborious.  With  the  incident  field  then  given 
by  (III. 40)  the  primary  field  on  the  aperture  is 

4>C  -  +  $2  +  **  +  4>2  +  0(k"3/2)  (IV. 11) 

G  d 

where  <j>  and  <j>  are  given  in  Equations  (11.46)  through  (11.53).  Plots  of 
the  amplitude  and  phase  shift  for  are  given  in  Figures  41  and  42. 
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V.  APPLICATION:  WAVEGUIDE  LENS  ANTENNAS 


From  the  onset  of  radio,  optical  techniques  have  been  applied  to  meet 
the  needs  of  microwave  communications.  The  simplest  and  most  fundamental 
of  these  problems  has  been  the  collimation  of  microwave  energy  from  a 
point  source.  It  was  readily  recognized  that,  as  in  optics,  the  solution 
was  a  standard  lens  built  of  dielectric  material.  These  lenses,  however. 


were  heavy  and  were  quickly  discarded  in  favor  of  reflectors.  But  in  1944, 

Kock  [18]  suggested  that  a  metal  guide  supporting  modes  can  be  used  as  a 

refracting  medium,  with  its  index  of  refraction  given  by  X/X  ,  where  X  is 

8  8 

the  guide  wavelength.  This  opened  the  possibility  of  significantly  reducing 
weight  and  cost,  and  because  of  the  less  stringent  manufacturing  tolerance, 
the  lens  began  to  compete  with  other  antenna  designs.  In  1950,  Ruze  [28] 


showed  that  a  "constrained"  lens,  that  is,  a  lens  operating  in  the  dominant 


TE1q  mode,  can  be  used  for  wide-angle  scanning  applications.  Such  a  two- 
dimensional  lens  is  shown  in  Figure  43.  The  guides  are  aligned  along  the 


contours  1^  and  Z^  with  the  edges  parallel  to  the  E  vector  so  that  only  the 
TE^q  mode  propagates.  Thus,  if  the  plate  surfaces  are  space  a  apart  we 
can  define  a  refractive  index  n  as 


(V.l) 


and  then  the  phase  shift  through  any  guide  of  length  t  is  27rnt/X.  Note 
that  in  contrast  to  natural  dielectrics  n  is  less  than  one.  In  analogy  to 
dielectric  lenses,  the  design  of  the  waveguide  lens  is  then  obtained  via 
standard  phase-path  length  analysis.  Figure  43  shows  that  for  a  feed  at 
a  focal  point  offset  a  degrees  from  the  z  axis  and  the  main  beam  in  the 
same  direction,  the  three  electrical  path  lengths  FPQN,  FOO’N',  and  FP'Q' 
are  set  equal  and  the  equations 
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z"  +  2zf  cos  a  +  cos^  a  -  0  (V.2) 

and 

nQtQ  +  nt  +  (t  -  tg  +  z)  cos  a  ■  0  (V.3) 

are  obtained.  Equation  (V.2)  gives  the  profile  of  the  first  surface  I.. 

This  is  an  ellipse  with  its  major  axis  perpendicular  to  the  lens  axis. 

Equation  (V.3)  is  one  relation  between  n  and  t  and  a  second  arbitrary 
relation  can  be  chosen.  Some  of  the  commonly  chosen  are 

(a)  a  lens  of  constant  thickness 

(b)  a  lens  with  a  plane  second  surface 

(c)  a  third  focal  point 

(d)  a  constant  refractive  index. 

There  are  two  main  features  that  distinguish  the  metal-plate  lens  from 
a  dielectric  lens.  First,  since  n  <  1,  E^  and  E 2  are  concave.  Second,  the 
rays  in  the  lens  are  "constrained."  By  this  it  is  meant  that  the  rays  in 
the  lens  always  travel  parallel  to  the  z  axis ,  regardless  of  the  angle  of 
incidence.  It  is  this  property  that  gives  these  lenses  such  good  wide-angle 
scanning  ability. 

The  standard  design  just  outlined  has  a  basic  limitation.  It  models 
the  region  of  any  waveguide  by  a  medium  with  index  of  refraction  n,  given 
by  (V.l),  and  does  not  consider  the  actual  physical  structure  and  the  scat¬ 
tering  off  the  plates.  It  assumes  that  the  phase  shift  at  E^  and  Z^  is  zero  or 
very  small,  but  more  fundamentally,  the  design  is  incomplete  as  it  cannot 
predict  the  amplitude  distribution  in  the  lens.  Note  that  because  the  lens 
is  "constrained,"  Snell's  Law  cannot  be  used  to  give  even  a  qualitative 
idea  of  what  the  amplitude  distribution  is.  What  is  typically  done,  as 
in  a  recent  design  by  Dion  and  Ricardi  [13],  is  to  assume  that  the  power 


1  - 
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transmission  coefficient  at  the  interface  of  any  guide  is  simply  one. 

With  the  methods  in  this  paper  the  amplitude  distribution  can  be  found. 

As  an  example,  we  study  the  scanner  in  Figure  44.  We  choose  2^  to  be 
a  plane,  and  ignoring  interactions  between  and  2^  and  any  interactions 
between  the  guides,  we  apply  the  methods  of  this  chapter  to  each  guide  in 
the  lens.  The  principles  of  operation  of  this  lens  are  well-known  [28],  and 
the  pertinent  parameters  are:  f  ■  14  X,  dg  ■  3.5  X,  n^  *  0.8,  a  ■  18°, 
f/D  *  0.95,  and  number  of  guides  *  21. 

Since  we  have  already  shown  that  the  phase  shift  at  is  minimal,  we 
concentrate  on  the  power  transmission.  Since  the  distance  from  the  feed  to 
any  waveguide  is  much  greater  than  the  width  of  the  guide,  we  approximate 
the  field  at  each  guide  by  a  plane  wave  and  examine  the  power  transmission 
coefficient  at  each  guide  as  the  feed  sweeps  along  the  focal  arc  shown  in 
Figure  44  .  The  use  of  a  plane  wave  allows  us  to  include  the  first  interaction 
when  computing  the  power  transmission  coefficients.  These  power  transmission 
coefficients  are  shown  in  Figures  45  thru  50  for  six  different  scan  angles. 

We  observe  that  for  small  scan  angles  the  amplitude  distribution  in  the  lens 
is  nearly  symmetrical  and,  except  near  the  edges,  has  an  additional  taper  of 
approximately  1  dB.  As  the  scan  angle  increases,  the  distribution  now  becomes 
asymmetrical.  The  staggers  in  the  lower  half  of  the  lens  channel  more 
power  into  the  guides  and  those  on  top  act  to  block  out  some  of  the  incident 
field.  For  scan  angles  approaching  a,  the  amplitude  distribution  becomes 
increasingly  flat  near  the  center  of  the  lens  with  variation  of  approximately 
2  dB  near  the  edges.  We  also  note  that  as  the  scan  angle  increases  less 
power  is  coupled  into  the  lens.  This  is  shown  in  Figure  51  where  we  plot 
the  efficiency  at  1^,  the  ratio  of  the  power  in  the  lens  to  the  incident 
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power,  as  a  function  of  scan  angle.  For  the  scan  angle  a  ■  18®,  the  efficiency 
is  down  by  0.7  dB  from  its  maximum  at  6  ■  0. 

The  effects  on  the  secondary  pattern  of  a  non-uniform  power  transmission 
thru  the  lens  are  examined  in  Figures  52  thru  58  .  We  compare  the  H-plane 
far  field  pattern  with  the  amplitude  distribution  calculated  by  the  methods 
in  this  thesis  with  the  pattern  given  by  assuming  a  power  transmission 
coefficient  of  one  through  each  guide,  as  done  by  Dion  and  Ricardi  [13].  Each 
pattern  is  normalized  by  the  maximum  intensity  computed  by  our  method.  We 
observe  that  the  nonuniform  power  transmission  distributions  in  Figures  52 
through  57  result  in  a  side-lobe  level  approximately  1  dB  below  that  calculated 
by  Dion  and  Ricardi.  Also,  the  intensity  of  the  main  beam  is  different. 

The  variation  of  the  main  beam  strength  with  the  scan  angle  is  given  in 
Figure  58.  The  curves  are  normalized  by  the  main  beam  strength  at  B  *  0, 
calculated  by  our  method.  We  note  that  the  main  beam  strength  has  a  stronger 
variation  than  would  be  predicted  by  Dion  and  Ricardi,  with  the  strength 
down  by  1  dB  at  B  -  14®. 
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Figure  56.  H-plane  pattern,  3  ■  16 


VI.  CONCLUSIONS 


In  chis  thesis,  we  examined  the  scattering  of  a  plane  and  cylindrical 
wave  by  a  staggered  parallel-plate  wave  guide,  shown  in  Figure  1.  In  partic¬ 
ular,  we  concentrated  on  the  fields  in  the  guide  for  0.5X<a<lA.  In 
Chapter  II,  a  quick  review  of  some  relevant  techniques  was  given,  and  in 
Chapter  III,  past  approaches  to  scattering  by  a  parallel-plate  guide  were  out¬ 
lined  and  critically  examined.  We  concluded  in  Chapter  III  that  the  Aperture 
Field  Method  offered  the  greatest  accuracy  and  flexibility  for  the  problem  at 
hand.  The  first  interaction  between  the  plates  was  also  exactly  formulated 
and  solved  for  in  Chapter  IV. 

In  Chapter  IV  we  used  the  Aperture  Field  Method  to  examine  the  TE^q 
mode  in  a  staggered  parallel-plate  guide.  We  gave  quantitative  results  for 
the  power  transmission  coefficient  and  showed  that  the  phase  shift  at  the 
aperture  of  the  guide  was  very  small  and  weakly  dependent  on  Because 

of  this  weak  dependence,  the  phase  shift  was  well-approximated  by  its  value 
at  1  *  0.  Tables  to  aid  in  this  computation  are  included  in  Appendix  B. 

In  Chapter  V,  we  applied  the  Aperture  Field  Method  to  the  study  of 
waveguide  lens  antennas  and  showed  how  the  amplitude  distribution  in  the 
lens  may  be  obtained,  something  previously  not  available. 

Finally,  we  noted  that,  although  we  have  concentrated  on  the  TE.^  mode 
and  the  E-wave,  the  same  methods  may  be  used  for  the  H-wave  problem.  Also, 
the  same  aperture  field  formulation  can  be  used  to  find  the  reflection 
coefficient  for  the  field  incident  from  inside  the  guide.  One  merely 
expresses  the  incident  field  as  a  sum  of  two  plane  waves  and  for  each 
plane  wave  makes  use  of  Equations  (IV. 3),  (IV. 6)  and  (IV. 7). 
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APPENDIX  A 


FRESNEL  INTEGRAL  WITH  COMPLEX  ARGUMENTS 
In  this  appendix,  some  basic  properties  of  the  Fresnel  integral 
with  complex  arguments  are  presented.  The  Fresnel  integral  is  defined 
as 


-iir/4  ®  2 

f«)  «•*-*-* -/„ 


(A.  1) 


where  the  integration  path  in  /  is  any  path  that  goes  to  infinity  in 

5 

the  first  quadrant.  The  asymptotic  expansion  of  F(£)  may  be  shown  as 


2  iir/4  00 

F(£)  -  0[-Re(£e"llr/4)]  +  e15  £  T (n  +  1/2) (i£2) 


-n 


n-0 


-  9(-Re(£e~ilT/4)]  +  F (£)  +  0(£-3  e15  )  for  |£|  »  0 


(A.  2) 


where  0  is  the  unit  step  function  and 


F(£) 


sir  +  iir/4 

2771 


(A.  3) 


Furthermore,  one  can  establish  the  following  inequality 


r,  if2 

7  I*15  I 


|f«)|  << 


vi  +  +  |.1( 


when  Re(£e  ilT^4)  >  0 


when  Re(£e  ^1T^4)  <  0 


(A. 4) 


which  can  be  derived  from  (Al)  by  introducing  the  change  of  variable 
t  ■  t  +  £.  Substitution  of  the  variable  -£  into  (Al)  allows  one  to 
arrive  at 

F(£)  +  F(-£)  -  1,  (A. 5) 


for  complex  £'s. 
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There  are  available  some  fast  routines  for  accurate  determination  of 
the  Fresnel  integral  (Al).  Application  of  one  of  these  routines  (see 
Table  A-l)  allows  one  to  evaluate  F(C)  for  real  or  complex  C's.  Some 
representative  results  are  shown  in  this  appendix.  Figure  A-l  is  the  plot 
of  the  magnitude  and  phase  of  F(C)  for  real  C's.  Figures  A-2,  A-3  and  A-4 
are  the  three-dimensional  plots  of  the  magnitude,  the  real  and  the  imaginary 
parts  of  F(C)  for  complex  C's  in  the  first  quadrant,  i.e.,  Re  E>o  and  Im  £>o. 
Figures  A-5,  A-6  and  A-7  are  the  counterparts  of  A-2,  A-4,  and  A-5, 
respectively,  for  complex  C's  in  the  fourth  quadrant,  i.e.,  Re  %>o  and 
Im  C^o.  Clearly,  because  of  (A5)  the  behavior  of  F(C)  in  the  third  and 
the  second  quadrants  is  similar  to  its  behavior  in  the  first  and  the  fourth 
quadrants,  respectively. 


TABLE  A-l 


FRESNEL  INTEGRAL  WITH  COMPLEX  ARGUMENT 


-iir/4 

ZFR-  V- 


dt 


SUBROUTINE  FRNL(Z1, ZFR) 

IMPLICIT  REAL* 8  (A-H,  O-Y) 

IMPLICIT  COMPLEX*16  (Z) 

REAL *8  IM 
PI-3.1415926535DO 
ZI*(ODO, IDO) 

SQPI-DSQRT(PI) 

Z-Z1*DSQRT (2 . DO/PI) 

ZW«SQPI*Z/ (l.-Zl) 

X-DREAL(ZW) 

Y-DIMAG(ZW) 

ABSX-DABS(X) 

ABSY-DABS(Y) 

C*  ZlST:  Z  IN  1ST  QUADRANT 
Z 1ST-DCMPLX (ABSX, ABSY) 

CALL  KOEZ (ABSX, ABSY, RE, IM) 

C*  1ST  QUADRANT:  W(Z) 

ZFR-DCMPLX(RE.IM) 

C*  3RD  QUADRANT:  W(-Z)»2*EXP(-Z**2)-W(Z) 
CALL  UNDERZ('OFF') 

IF( (X.LE.O. ) .AND. (Y. LE.O. ) ) 
1ZFR“2.*CDEXP(-Z1ST*Z1ST) -ZFR 
C*  4TH  QUADRANT:  W(CONJ(Z))»CONJ(W(-Z)) 
IF((X.GE.O.) .AND. (Y.LE.O.)) 

IZFR-DCON JG ( 2 . *CDEXP ( -Z1ST*Z 1ST) -ZFR) 
C*  2ND  QUADRANT:  W(CONJ(-Z))-CONJ(W(Z)) 
IF( (X.LE.O.). AND. (Y.GE.O.)) 
lZFR-DCONJG(ZFR) 

ZFR- ( 1 . +ZI) / 2 . *( 1 . -ZFR*CDEXP ( ZW*ZW) ) 
CALL  UNDERZ ( ' ON ' ) 

ZFR- . 5DO- ( 1 . DO-ZI ) *ZFR/ 2 . DO 

RETURN 

END 


SUBROUTINE  WOFZ(X,  Y,  RE ,  IM) 

R£AL*8  X,Y,RE,IM 

INTEGER*4  CAPN,NU,N,NPI 

REAL *8  H, H2, LAMBDA, R1,R2,S, SI, S2,T1,T2,C 

LOGICAL* 4  B 

IF( (Y .LT . 4. 29) .AND . (X. LT. 5 . 33) )  GO  TO  1 
H-O. 

CAPN-0 
NU-8 
GOTO  2 

1  CONTINUE 

S-(1.-Y/4.29)*DSQRT(1.-X*X/28.41) 

H=  1 . 6*S 

H2-2.*H 

CAPN-6+23*S 

LAMBDA-H2**CAPN 

NU-9+21*S 

2  CONTINUE 

C*  IN  THE  FOLLOWING  STATEMENT,  LAMB DA *0  COVERS  THE  UNDERFLOW  CASE  WHERE 
C*  H.GT.O.  IS  VERY  SMALL. 

B-((H.EQ.O.) .OR. (LAMBDA. EQ.O.)) 

Rl-O. 

R2-0 . 

51- O. 

S2=0. 

N-NU 

3  CONTINUE 
NP1-N+1 
T1=Y+H+NP1*R1 
T2-X-NP1*R2 
C-0.5/(Tl*Tl+T2*T2) 

Rl-C*T1 

R2-C*T2 

IF( (H.GT.O.) .AND. (N.LE.CAPN))  GOTO  4 
GOTO  5 

4  CONTINUE 
T1-LAMBDA+S1 

S  1-R1*T1-R2.*S  2 

52- R2*Tl+Rl*S2 
LAMBDA-LAMBDA/H2 

5  CONTINUE 
N-N-l 

IF(N.GE.O)  GOTO  3 
IF(Y.EQ.O.)  RE-DEXP(-X*X) 

IF((Y.NE.O.) .AND.B)  RE-  1 . 12837916709551*R1 

IF( (Y.NE.O. ) .AND. ( .NOT ,B) ) RE-  1.12837916709551*S1 

IF(B)  IM  -  1 . 1 283 791 6709 551*R2 

IF ( . NOT. B)  IM-  1 . 12837916709551*S2 

RETURN 
END 


Figure  A-6.  Real  part  of  F(5),  defined  in  (Al),  for  £'s  in  the  fourth  quadrant 
of  complex  £-plane.  Re  F(5)  is  truncated  at  -5  or  5. 


Figure  A-7.  Imaginary  part  of  F(?),  defined  in  (Al) ,  for  £'s  in  the  fourth 
quadrant  of  complex  5-plane.  Im  F(5)  is  truncated  at  -5  or  5. 


APPENDIX  B 


THE  FUNCTION  G  (a) 

As  illustrated  in  this  report,  a  central  step  in  the  Wiener-Hopf 
Method  as  applied  to  the  waveguide  is  the  computation  of  the  "plus" 
part  of 


G(a) 


1  +  e 


-2(a2-k2)1/2b 


(B.l) 


where  for  the  latter  case  (+) ,  we  say  G(a)  *  G(a).  One  expression  for 
G+.(a) ,  already  presented  in  Equations  (III. 14b)  and  (III. 14c), 
is  the  infinite  product  form.  For  small  or  moderate  guide  dimensions, 
this  product  converges  quite  rapidly;  however,  for  larger  guide  dimensions, 
the  infinite  product  form  can  be  very  laborious.  As  a  matter  of  fact,  it 
can  be  shown  that  for  truncation  of  the  product  after  N  terms, the  error, 

E,  is 

E  -  [  (ab/ir)  2/N]  if  e  «  1  .  (B.2) 


Then  for  ab  -  not  small  -  N  must  be  exceedingly  large. 
Another  useful  formula  for  G+(a)  is 


G+(a)  »  exp[U(a/k,  2kb)  ]  a  >  0^"  , 


(B.3) 


where  U(s,p)  is  the  Weinstein  function  [21],  [31],  given  by 

OCs.p)  -  | 


-  2  212  1/2 
“  ln(l  +  exp[ip  -  pt4]) (1  +  it  ) (1  +  it  ;  dt 


t(l  +  \  it2)  -  2  ^2  s  exp[i7r/4] 


(B.4) 


Equations  (B.7)  and  (B.8)  in  [21]  are  only  valid  for  a  >  0. 


>  iu'| 


T 
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the  plus  in  ln(l  +  exp[*])  being  used  for  G(a)  -  G(a). 
For  a  <  0,  this  relationship  is  used: 


G+(a)  -  (e~ikb  +  elkb) l/2eikb/2 

exp  2^1  {  m{l  + 

P 

where  the  path  P  goes  from  -®>  to  0+  below  the  real  axis  of  s,  circles 
around  the  origin,  and  then  goes  from  0+  to  +«  above  the  real  axis. 

Note  that  the  integral  in  (B.7)  converges  quite  rapidly.  Other  asymptotic 
formulas  for  special  cases  may  be  found  in  the  references  listed  above. 

Tables  B.l  and  B.2  and  Figures  B.l  and  B.2  give  the  values  of  G+(a) 
for  different  values  of  kb. 


J 


114 


TABLE  B.l 

Values  of  G+(kx),  (Magnitude,  Phase). 

....  2  ..1/2 

G(kx)  -  1  -  e"2kb(x  _1) 


Values  of  G+(kx)  for  kb  *»  1.0 

Values  of  G,(kx)  for  kb  —  2.0 

X 

G.(kx) 

X 

G,(kx) 

X 

G.(kx) 

X 

GJkx) 

-1.00 

(0.00, 

0.0) 

0.00 

(1.30, 

-16.4) 

-1.00 

(0.00, 

0.0) 

0.00 

(1.35, 

12.3) 

-0.95 

(0.50, 

-67.6) 

0.05 

(1.30, 

-15.3) 

-0.95 

(1.12, 

-63.1) 

0.05 

(1-32, 

12.6) 

-0.90 

(0.69, 

-60.2) 

0.10 

(1.30, 

-14.3) 

-0.90 

(1.46, 

-49.3) 

0.10 

(129, 

12.9) 

-0.85 

(0.82, 

-54.7) 

0.15 

(1.29, 

-13.4) 

-0.85 

(1.64, 

-39.2) 

0.15 

(1.26, 

13.0) 

-0.80 

(092, 

-50.1) 

0.20 

(1.29, 

-12.5) 

-0.80 

(1.75, 

-31.1) 

0.20 

(1.23, 

13.0) 

-0.75 

(0.99, 

-46.2) 

0.25 

(1.29, 

-11.7) 

-0.75 

(1.81, 

-24.4) 

0.25 

(1.21, 

12.9) 

-0.70 

(1.05, 

-42.8) 

0.30 

(1.28, 

-10.9) 

-0.70 

(1.84, 

-18.7) 

0.30 

(1.19. 

12.8) 

-0.65 

(1.10, 

-39.7) 

0.35 

(1.28, 

-10.2) 

-0.65 

(1.84, 

-13.8) 

0.35 

(1-17. 

12.6) 

-0.60 

(1.14, 

-36.9) 

0.40 

(1.27, 

-9.6) 

-0.60 

(1.82, 

-9.5) 

0.40 

(1.15, 

12.4) 

-0.55 

(1.18, 

-34.4) 

0.45 

(1.27, 

-8.9) 

-0.55 

(1.79, 

-5.8) 

0.45 

(114, 

12.1) 

-0.50 

(1.20, 

-32.0) 

0.50 

(1.27, 

-8.3) 

-0.50 

(1.76, 

-2.6) 

0.50 

(1.12, 

11.8) 

-0.45 

(1.23, 

-29.9) 

0.55 

(1.26, 

-7.8) 

-0.45 

(1.72, 

0.2) 

0.55 

(1.11. 

11.5) 

-0.40 

(1.24, 

-27.9) 

0.60 

(1.25, 

-7.3) 

-0.40 

(1.68, 

2.6) 

0.60 

(1-10, 

11.2) 

-0.35 

(1.26, 

-26.1) 

0.65 

(1.25, 

-6.8) 

-0.35 

(1.63, 

4.7) 

0.65 

(1-09, 

10.9) 

-0.30 

(1.27, 

-24.4) 

0.70 

(1.24, 

-6.3) 

-0.30 

(1.59, 

6.5) 

0.70 

(1.08, 

10.5) 

-0.25 

(1-28, 

-22.8) 

0.75 

(1.24, 

-5.9) 

-0.25 

(1.54, 

8.0) 

0.75 

(1.07, 

10.2) 

-0.20 

(129, 

-21.3) 

0.80 

(1-23, 

-5.5) 

-0.20 

(1.50, 

9.3) 

0.80 

(1-06, 

9.9) 

-0.15 

(1.29, 

-20.0) 

0.85 

(1.23, 

-5.1) 

-0.15 

(1.46, 

10.3) 

0.85 

(1.06, 

9.5) 

-0.10 

(1.29, 

-18.7) 

0.90 

(1.22, 

-4.8) 

-0.10 

(1.42, 

11.2) 

0.90 

(1.05, 

9.2) 

-0.05 

(1.30, 

-17.5) 

0.95 

(1.22, 

-4.5) 

-0.05 

(1.38, 

11.8) 

0.95 

(1.04, 

8.9) 

1.00 

(1.21, 

-4.2) 

1.00 

(1.04, 

8.6) 

Values  of  G+(kx)  for  kb  ••  3.0  Values  of  Gr(kx)  for  kb  **  3.1416 


X 

G.(kx) 

X 

G.(kx) 

X 

GAkx) 

X 

G.(kx) 

-1.00 

(0.00, 

0.0) 

0.00 

(0.53, 

40.9) 

-1.00 

(0.00, 

0.0) 

0.00 

(0.00, 

135.0) 

-0.95 

(2.11. 

-35.2) 

0.05 

(0.51, 

34.6) 

-0.95 

(1.96, 

162.7) 

0.05 

(0.08, 

137.7) 

-0.90 

(2.56, 

-14.1) 

0.10 

(0.51, 

28.5) 

-0.90 

(2.35, 

-174.3) 

0.10 

(0.14, 

140.3) 

-0.85 

(2.70, 

1.3) 

0.15 

(0.51, 

23.0) 

-0.85 

(2.44, 

-157.2) 

0  15 

(0.20, 

142.6) 

-0.80 

(2.67, 

13.7) 

0.20 

(0.52, 

18.2) 

-0.80 

(2.39, 

-143.2) 

0.20 

(025, 

144.8) 

-0.75 

(2.56, 

23.9) 

0.25 

(0.53, 

14.2) 

-0.75 

(2.25, 

-131.3) 

0.25 

(0.30. 

146.8) 

-0.70 

(2.40, 

32.5) 

0.30 

(0.55, 

11.0) 

-0.70 

(2.08, 

-120.8) 

0.30 

(0.34, 

148.7) 

-0.65 

(2.22. 

39.9) 

0.35 

(0.57, 

8.4) 

-0.65 

(I. 88, 

-111.6) 

0.35 

(0.38, 

150.4) 

-0.60 

(2.02, 

46.1) 

0.40 

(0.59, 

6.3) 

-0.60 

(1.68, 

-103.3) 

0.40 

(0.41, 

152.0) 

-0.55 

(1.82, 

51.3) 

0.45 

(0.61, 

4.6) 

-0.55 

(1.48, 

-95.8) 

0.45 

(0.44, 

153.5) 

-0.50 

(1.64, 

55.6) 

0.50 

(0.63, 

3.3) 

-0.50 

(1.28, 

-88.9) 

0.50 

(0.47, 

154.8) 

-0.45 

(1.46, 

58.9) 

0.55 

(0.65, 

2.2) 

-0.45 

(1.10, 

-82.7) 

0.55 

(0.49, 

156.1) 

-0.40 

(1.29, 

61.3) 

0.60- 

•  (0.67, 

1.4) 

-0.40 

(0.93, 

-77.0) 

0.60 

(0.51, 

157.3) 

-0.35 

(1.14, 

62.7) 

0.65 

(0.68, 

0.7) 

-0.35 

(0.77, 

-71.S) 

0.65 

(0.53, 

15S.4) 

-0.30 

(1-00, 

63.0) 

0.70 

(0.70, 

0.2) 

-0.30 

(0.63, 

-67.0) 

0.70 

(0.55, 

159.4) 

-0.25 

(0.88, 

62.2) 

0.75 

(0.71, 

-0.2) 

-0.25 

(0.49, 

-62.5) 

0.75 

(0.57, 

160.3) 

-0.20 

(0.77, 

60.2) 

0.80 

(0.73, 

-0.5) 

-0.20 

(0.37, 

-58.4) 

0.80 

(0.5S. 

161.2) 

-0.15 

(0.69, 

57.0) 

0.85 

(0.74, 

-0.8) 

-0.15 

(0.26, 

-54.7) 

0.85 

(0.60, 

162.0) 

-0.10 

(0.62, 

52.5) 

0.90 

(0.75, 

-1.0) 

-0.10 

(0.16, 

-51.2) 

0.90 

(0.61, 

162.8) 

-0.05 

(0.57, 

47.1) 

0.95 

(0.76, 

-1.1) 

-0.05 

(0.08, 

-4S.0) 

0.95 

(0.62. 

163  5) 

1.00 

(0.77, 

-1-3) 

1.00 

(063, 

164.1) 

-  -  jggr 
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TABLE  B.l  (Continued) 


Values  of  Gr(kx)  for  kb  — 

3.2 

Values  of  G.(kx)  for  kb  — 

4.0 

X 

G.(kx) 

X 

G,(kx) 

X 

GJkx) 

X 

G.i.kx) 

-1.00 

(0.00, 

0.0) 

0.00 

(0.34, 

136.7) 

-1.00 

(0.00, 

0.0) 

0.00 

(123, 

159.6) 

-0.95 

(1.93, 

162.7) 

0.05 

(0.41, 

139.4) 

-0.95 

(1.69, 

164.7) 

0.05 

(1.23, 

161.9) 

-0.90 

(2.29, 

-173.9) 

0.10 

(0.47, 

141.9) 

-0.90 

(1.75, 

-166.7) 

0.10 

(1.24, 

163.9) 

-0.85 

(2.33. 

-156.6) 

0.15 

(0.52, 

144.2) 

-0.85 

(1.52, 

-145.6) 

0.15 

(1.23, 

165.6) 

-0.8  0 

(2.23, 

-142.4) 

0.20 

(0.57, 

146.4) 

-0.80 

(1.19, 

-128.6) 

020 

(1.23, 

167.2) 

-0.75 

(2.06, 

-130.2) 

0.25 

(061, 

148.4) 

-0.75 

(0.83, 

-114.2) 

0.25 

(1.22, 

168.5) 

-0.70 

(1.85, 

-119.7) 

0.30 

(0.64, 

150.2) 

-0.70 

(0.49, 

-101.7) 

0.30 

(1.21. 

169  8) 

-0.65 

(1.62, 

-110.3) 

0.35 

(0.68, 

151.9) 

-0.65 

(0.18, 

-90.8) 

0.35 

(1.21, 

170.8) 

-0.60 

(1.39, 

-101.9) 

0.40 

(0.70, 

153.4) 

-0.60 

(0.10, 

98.9) 

0.40 

(120, 

171.7) 

-0.55 

017, 

-94.3) 

0.45 

(0.73, 

154.9) 

—0.55 

(0.34, 

107.5) 

0.45 

(1.19, 

172.6) 

-0.50 

(0.96, 

-87.4) 

0.50 

(0.75, 

156.2) 

-0.50 

(0.54, 

iu.2) 

0.50 

(U8, 

173.3) 

-0.45 

(0.77, 

-81.2) 

0.55 

(0.77, 

157.4) 

-0.45 

(0.70, 

122.1) 

0.55 

(1.17, 

173.9) 

-0.40 

(0.59. 

-75.4) 

0.60 

(0.79, 

158.6) 

-0.40 

(0.84, 

128.3) 

0.60 

(1.17, 

174.5) 

-0.35 

(0.42. 

-70.1) 

0.65 

(0.80, 

159.6) 

-0.35 

(0.94, 

133.9) 

0.65 

(1.16. 

174.9) 

-0.30 

(0.28, 

-65.3) 

0.70 

(0.82, 

160.6) 

-0.30 

(103, 

138.9) 

0.70 

(115. 

175.4) 

-0.25 

(0.14, 

-60.8) 

0.75 

(0.83, 

161.5) 

-0.25 

(1.09, 

143.4) 

0.75 

(1.14, 

175.7) 

-0.20 

(0.02. 

-56.8) 

0.80 

(0.84, 

162.4) 

-0.20 

(1.14, 

147.4) 

0.80 

(1.14. 

176.1) 

-0.15 

(0.09, 

127.0) 

0.85 

(0.85, 

163.1) 

-0.15 

(1.18, 

151.0) 

0.85 

(1.13. 

176.4) 

-0.10 

(0.18, 

130.5) 

0.90 

(0.86, 

163.9) 

-0.10 

(1.20, 

154.2) 

0.90 

(1.12, 

176.6) 

-0.05 

(0.27, 

133.7) 

0.95 

(0.87, 

164.5) 

-0.05 

(1.22, 

157.0) 

0.95 

(1.12, 

176.9) 

1.00 

(0.87, 

165.2) 

1.00 

(111. 

177.1) 

Values  of  G.(kx)  for  kb  =•  5.0 

Values  of  Gt(kx)  for  kb  »• 

6.0 

X 

G.(kx) 

X 

G.(kx) 

X 

G.(kx) 

X 

G,(kx) 

-1.00 

(0.00, 

0.0) 

0.00 

(1.38,  -171.8) 

-1.00 

(0.00, 

0.0) 

0.00 

(0.75,  -143.1) 

-0.95 

(1.92, 

174.0) 

0.05 

(1.34,  -171.0) 

-0.95 

(2.22. 

-164.6) 

0.05 

(0.71,  -148.1) 

-0.90 

(1.57, 

-150.8) 

0.10 

(1.30,  -170.5) 

-0.90 

(1.18, 

-122.3) 

0.10 

(0  69,  -153.1) 

-0.85 

(0.93, 

-125.0) 

0.15 

(1.26,  -170.2) 

-0.85 

(0.05, 

88.7) 

0.15 

(0.68,  -157.7) 

-0.80 

(0.27, 

-104.3) 

0.20 

(1.22,  -170.2) 

-0.80 

(106, 

113.6) 

0.20 

(0  68,  -161.7) 

-0.75 

(0.31, 

93.0) 

0.25 

(1.19,  -170.3) 

-0.75 

(1.78. 

134.4) 

0.25 

(0.69,  -165.0) 

-0.70 

(0.78. 

107.8) 

0.30 

(1.17,  -170.4) 

-0.70 

(2.23, 

152.2) 

0  30 

(0.71,  -167.7) 

-0.65 

(1.14, 

120.5) 

0.35 

(1.15,  -170.7) 

-0.65 

(2.46, 

167.4) 

0.35 

(0.72,  -169.8) 

-0.60 

(1.40, 

131.7) 

0.40 

(1.13,  -171.0) 

-0.60 

(2.51, 

-179.4) 

0.40 

(0.74,  —171.5) 

-0.55 

(1.58, 

141.4) 

0.45 

(1.11,  -171.4) 

-0.55 

(2.45, 

-168.0) 

0.45 

(0.75,  -172.9) 

-0.50 

(1.68, 

149.9) 

0.50 

(1.10,  -171.8) 

-0.50 

(2.31, 

-158.3) 

0.50 

(0.77,  -174.0) 

-0.45 

(1.74, 

157-2) 

0.55 

(1.09,  -172.1) 

-0.45 

(2.12, 

-150.1) 

0.55 

(0.78,  -174.9) 

-0.40 

(1.75, 

163.6) 

0.60 

(1.08,  -172.5) 

-0.40 

(1.92, 

-143.4) 

0.60 

(0.79,  —175.6) 

-0.35 

(1.74, 

169.1) 

0.65 

(1.07,  -172.8) 

-0.35 

(171, 

-138.1) 

0.65 

(0.80,  -176.2) 

-0.30 

(1.71, 

173.7) 

0.70 

(1.06,  -173.2) 

-0.30 

(1.50, 

-134.4) 

0.70 

(0.81,  -176.7) 

-0.25 

(1.66, 

177.6) 

0.75 

(1.06,  -173.5) 

-0.25 

(1.32, 

-132.2) 

0.75 

(0.82,  -177.0) 

-0.20 

(1.61, 

-179.1) 

0.80 

(1.05,  -173.8) 

-0.20 

(1.15, 

-131.5) 

0.80 

(0  83,  -177.4) 

—0.15 

(1.55, 

-176.5) 

0.85 

(1.05,  -174.1) 

-0.15 

(1.01, 

-132.4) 

0.85 

(0.84,  -1776) 

-0.10 

(1.49, 

-174.5) 

0.90 

(1.04,  -174.3) 

-0.10 

(0.90, 

-134.8) 

0.90 

(0.85,  -177  8) 

-0.05 

(1.44, 

-172.9) 

0.95 

(1.04,  -174.6) 

-0.05 

(0.81, 

-138.5) 

0.95 

(0.85,  -178.0) 

1.00 

(1.04,  -174.8) 

1.00 

(0.86,  -178.2) 
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TABLE  B.l  (Continued) 


Values  of  G*(fce)  for  kb  -  6.2832 

Values  of  G+(kx)  for  kb  —  7.0 

X 

G.(kx) 

X 

GAkx) 

X 

Gt(kx) 

X 

G.(kx) 

-1.00 

(0.00, 

0.0) 

0.00 

(0.01, 

135.0) 

-1.00 

(0.00, 

0.0) 

0.00 

(115. 

155.5) 

-0.95 

(1.97, 

-145.5) 

0.05 

(0.10, 

138.9) 

-0.95 

(1.50, 

-141.2) 

0.05 

(1.16, 

158.7) 

-0.90 

(085. 

-100.4) 

0.10 

(0.19, 

142.4) 

-0.90 

(0.16, 

-91.4) 

0.10 

(1.17, 

161.4) 

-0.85 

(0.37, 

113.0) 

0.15 

(0.26, 

145.5) 

-0.85 

(0.95, 

125.4) 

0.15 

(1.18, 

163.7) 

-0.80 

(1.31, 

140.1) 

0.20 

(0.33. 

148.3) 

-0.80 

(1.58, 

155.0) 

0.20 

(1.17, 

165.7) 

-0.75 

(1.93, 

163.0) 

0.25 

(0.38, 

150.3) 

-0.75 

(1.80, 

180.0) 

0.25 

(1.17, 

167.4) 

-0.70 

(2.26, 

-177.2) 

0.30 

(0.43, 

153.0) 

-0.70 

(1.72, 

-158.5) 

0.30 

(116, 

168.9) 

-0.65 

(2.38, 

-159.8) 

0.35 

■  (0.47, 

155.0) 

-0.65 

(1.47, 

-139.7) 

0.35 

(1.16, 

170.1) 

-0.60 

(2.33, 

-144.3) 

0.40 

(0.50, 

156.8) 

-0.60 

(1.12, 

-123.2) 

0.40 

(1.15, 

171.1) 

-0.55 

(2.18, 

-130.5) 

0.45 

(0.53, 

158.4) 

-0.55 

(0.75, 

-108.5) 

0.45 

(1.14, 

172.0) 

-0.50 

(1.97, 

-118.1) 

0.50 

(0.56, 

159.3) 

-0.50 

(0.39, 

-95.5) 

0.50 

(1.13, 

172.8) 

-0.45 

(1.72, 

-106.9) 

0.55 

(0.58, 

161.1) 

-0.45 

(0.06, 

-83.9) 

0.55 

(1.13, 

173.5) 

-0.40 

(1.47, 

-96.8) 

0.60 

(0.60, 

162.3) 

-0.40 

(0.23, 

106.4) 

0.60 

(1.12, 

174.0) 

-0.35 

(1.22, 

-87.8) 

0.65 

(0.61. 

163.3) 

-0.35 

(0.47, 

115.6) 

0.65 

(1.11. 

174.5) 

-0.30 

(0.98, 

-'9.6) 

0.70 

(0.63, 

164.3) 

-0.30 

(0.66, 

123.7) 

0.70 

(1.11, 

174.9) 

-0.25 

(0.77, 

-72. 2) 

0.75 

(0.64, 

165.1) 

-0.25 

(0.81, 

130.9) 

0.75 

(110, 

175.3) 

-0.20 

(0.57, 

-65.6) 

0.80 

(0.66, 

165.9) 

-0.20 

(0.93. 

137.2) 

0.80 

(1.10 

175.6) 

-0.15 

(0.40, 

-59.6) 

0.85 

(0.67, 

166.6) 

-0.15 

(101, 

142.8) 

0.85 

(1.09, 

175.9) 

-0.10 

(0.24, 

-54.2) 

0.90 

(0.68, 

167.3) 

-0.10 

(1.07, 

147.7) 

0.90 

(1.09, 

176.2) 

-0.05 

(0.11. 

-49.3) 

0.95 

(0.68, 

167.9) 

-0.05 

(1.12, 

151.9) 

0.95 

(1.08, 

176.4) 

1.00 

(0.69, 

168.4) 

1.00 

(1.08, 

176.6) 

TABLE  B . 2 


Values  of  G+(kx) ,  (Magnitude,  Phase) 
„/,_x  _  ,  ,  -2kb(x2-i)17 


V»lu«t  of  G+(lue)  for  leb  - 


12.33, 

-12.75) 

(2.22, 

4.53) 

(2.11, 

10.98) 

(2.01, 

15.52) 

(1.92, 

19.02) 

(1.83, 

21.81) 

(1.74, 

24.07) 

(1.67. 

25.90) 

(1.59, 

27.38) 

(1.53. 

28.56) 

(1.46. 

29.46) 

(1.40, 

30.13) 

(1.35. 

30.581 

(1.30, 

30.84) 

(1.25. 

30.92) 

(1.21, 

30.351 

(1.17. 

30.63) 

(1.13, 

30.291 

(1.10. 

29.34) 

(1.07, 

29.28) 

-168.47) 

-132.18) 

-117.47) 

-106.48) 

-97.47) 

-89.76) 

-33.01) 

-77.00) 

108.40) 

113.31) 

117.78) 

121.37) 

125.66) 

129.14) 

132.35) 

135.33) 

138.09) 

140.65) 

143.03) 

145.24) 


(0.36,  13.36) 


(0. 

91, 

147, 

,30) 

(0. 

.94, 

149. 

.21) 

(0. 

.96, 

150. 

,99) 

(0. 

,98, 

152. 

,64) 

(1. 

,00, 

154, 

.19) 

(1. 

,01, 

155, 

.62) 

(1. 

,02, 

156, 

.96) 

(1. 

,03, 

158, 

.21) 

(1. 

,04, 

159, 

.37) 

(1. 

,05, 

160, 

,45) 

(1. 

,05, 

161, 

.45) 

(1. 

,06, 

162, 

.39) 

(1. 

,06, 

163, 

.27) 

(1. 

,07, 

164, 

.09) 

(1. 

,07, 

164, 

.35) 

(1. 

,07, 

165 , 

,56) 

(1. 

.07, 

166, 

.22) 

(1. 

.07, 

166, 

.34) 

(1. 

,07, 

167, 

.  -2) 

(1. 

.07. 

167 

.96) 

(1.07,  168.47) 

TABLE  B.2  (Continued) 


Valuta 

of  G^tkx) 

for  Scb  “  3.0 

X 

G+(Vx) 

X 

«♦ 

(Sex) 

-1.00 

(1.79, 

176.94) 

0.00 

(1.41, 

175.94) 

-0.95 

(1.06, 

-129.43) 

0.05 

(1.38, 

177 . 26) 

-0.90 

(0.46, 

-108.21) 

0.10 

(1.36, 

178.39) 

-0.35 

(0.02, 

87.39) 

0.15 

(1.34, 

179.34) 

-0.80 

(0.40, 

99.96) 

0.20 

(1.32, 

-179.36) 

-0.75 

(0.70, 

110.53) 

0.25 

(1.30, 

-179.19) 

-0.70 

(0.93, 

119.62) 

0.30 

(1.28, 

-178.53) 

-0.65 

(1.11. 

127.55) 

0.35 

(1.26, 

-178.18) 

-0.60 

(1.24, 

134.52) 

0.40 

(1.24, 

-177.32) 

-0.55 

(1.34, 

140.69) 

0.45 

(1.23, 

-17'. 53) 

-0.50 

(1.41, 

146.16) 

0.50 

(1.21. 

-177.31) 

-0.45 

(1.46, 

151.03) 

0.55 

(1.20, 

-177.14) 

-0.40 

(1.49, 

155.36) 

0.60 

(1.19, 

-177.01) 

-0.35 

(1.50, 

159.20) 

0.65 

(1.17, 

-176.92) 

-0.30 

(1.51, 

162.60) 

0.70 

(1.16. 

-176.37) 

-0.25 

(1.50. 

165.62) 

0.75 

(1.15. 

-176.34) 

-0.20 

(1.49, 

168.27) 

0.30 

(1.14, 

-176.33) 

-0.15 

(1.47, 

170.60) 

0.35 

(1.14, 

-176.34) 

-0.10 

(1.45, 

172.64) 

0.90 

(1.13. 

-176.36) 

-0.05 

(1.43, 

174.41) 

0.95 

(1.12, 

-176.90) 

1.00 

(1.12. 

-176.94) 

Valuas 

of  G+(tac) 

for  Seb  “  3.1416 

X 

G^OiUC) 

X 

(Sue) 

-1.00 

(1.32, 

175.70) 

0.00 

(1.41, 

-130.00) 

-0.95 

(1.01, 

-128.21) 

0.05 

(1.39, 

-178.35) 

-0.90 

(0.36, 

-106.06) 

0.10 

(1.36, 

-177.91) 

-0.85 

(0.15, 

90.19) 

0.15 

(1.33, 

-177.14) 

-0.80 

(0.55, 

103.27) 

0.20 

(1.31, 

-176.52) 

-0.75 

(0.86, 

114.23) 

0.25 

(1.28, 

-176.04) 

-0.70 

(1.09, 

123.64) 

0.30 

(1.26, 

-175.67) 

-0.65 

(1.27, 

131.82) 

0.35 

(1.24, 

-175.39) 

-0.60 

(1.39, 

138.99) 

0.40 

a. 22, 

-175.20) 

-0.55 

(1.48, 

145.31) 

0.45 

(1.21, 

-175.07) 

-0.50 

(1.53, 

150.39) 

0.50 

(1.19, 

-175.00) 

-0.45 

(1.57, 

155.82) 

0.55 

(1.13, 

-174.98) 

-0.40 

(1.58, 

160.17) 

0.60 

(1.16, 

-174.49) 

-0.35 

(1.58, 

164.01) 

0.65 

(1.15, 

-175.03) 

-0.30 

(1.57, 

167.38) 

0.70 

(1.14, 

-173.09) 

-0.25 

(1.55, 

170.32) 

0.75 

(1.13, 

-175.17) 

-0.20 

(1.53, 

172.39) 

0.30 

(1.12, 

-175.27) 

-0.13 

(1.50, 

175.10) 

0.35 

(1.12, 

-175.37) 

-0.10 

(1.47, 

177.01) 

0.90 

(1.11, 

-175.43) 

-0.05 

(1.44, 

178.63) 

0.95 

(1.10, 

-175.59) 

1.00 

(1.10. 

-175.70) 

TABLE  B.2  (Continued) 


V»lu«s  of  G+(Scx)  for  lob  -  3.2000 

„ 

X 

G+(lcx)  * 

G+(Scx) 

-1.00 

(1.84, 

175.25) 

0.00 

(1.41, 

-178.33) 

-0.95 

(0.99, 

-127.64) 

0.05 

(1.38, 

-177.26) 

-0.90 

(0.32, 

-105.11) 

0.10 

(1.35, 

-176.40) 

-0.85 

(0.21, 

91.41) 

0.15 

(1.32, 

-175.71) 

-0.30 

(0.61, 

104.70) 

0.20 

(1.30, 

-175.18) 

-0.75 

(0.93, 

115.83) 

0.25 

(1.27, 

-174.73) 

-0.70 

(1.16. 

125.36) 

0.30 

(1.25. 

-174.49) 

-0.65 

(1.33, 

133.65) 

0.35 

(1.23, 

-174.29) 

-0.60 

(1.45, 

140.90) 

0.40 

(1.21, 

-174.18) 

-0.55 

(1.53, 

147.28) 

0.45 

(1.19, 

-174.12) 

-0.50 

(1.58, 

152.90) 

0.50 

(1.18, 

-174.12) 

-0.45 

(1.61, 

157.85) 

0.55 

(1.17, 

-174.15) 

-0.40 

(1.62, 

162.22) 

0.60 

(1.15, 

-174.22) 

-0.35 

(1.61, 

166.04) 

0.65 

(1.14, 

-174.32) 

-0.30 

(1.59, 

169.39) 

0.70 

(1.13, 

-174.43) 

-0.25 

(1.57, 

172.30) 

0.75 

(1.12, 

-174.55) 

-0.20 

(1.54, 

174.82) 

0.80 

(1.11. 

-174.69) 

-0.15 

(1.51, 

176.99) 

0.35 

(1.11, 

-174.33) 

-0.10 

(1.48, 

178.33) 

0.90 

(1.10, 

-174.97) 

-0.05 

(1.45, 

-179.62) 

0.95 

(1.09, 

-175.11) 

1.00 

(1.C9, 

-175.25) 

V»lu«»  of  G+(!tx)  for  Seb  »  4.0 


X 

K 

(lot) 

X 

i 

(Sex) 

-1.00 

(2.14, 

173.44) 

0.00 

(1.14, 

-155.41) 

-0.95 

(0.68, 

-115.37) 

0.05 

(1.10, 

-156.33) 

-0.90 

(0.37, 

92.56) 

0.10 

(1.06, 

-157.46) 

-0.35 

(1.10. 

112.93) 

0.15 

(1.02, 

-158.73) 

-0.80 

(1.59, 

129.21) 

0.20 

(1.00, 

-160.06) 

-0.75 

(1.89, 

142.73) 

0.25 

(0.98, 

-161.41) 

-0.70 

(2.07, 

154.19) 

0.30 

(0.96, 

-162.73) 

-0.65 

(2.15, 

164.00) 

0.35 

(0.95, 

-163.99) 

-0.60 

(2.15, 

172.42) 

0.40 

(0.94, 

-165.18) 

-0.55 

(2.11, 

179.64) 

0.43 

(0.94, 

-166.29) 

-0.50 

(2.04, 

-174.22) 

0.50 

(0.93, 

-167.31) 

-0.45 

(1.94, 

-169.04) 

0.55 

(0.93, 

-163.23) 

-0.40 

(1.84, 

-164.  77) 

0.60 

(0.93. 

-169.07) 

-0.35 

(1.73, 

-161.32) 

0.65 

(0.93, 

-169.84) 

-0.30 

(1.62, 

-158.65) 

0.70 

(0.93, 

-170.52) 

-0.25 

(1.52, 

-156.69) 

0.75 

(0.93, 

-171.14) 

-0.20 

(1.43, 

-155.39) 

0.30 

(0.93, 

-171. '0) 

-0.15 

(1.34, 

-154.68) 

0.35 

(0.93, 

-172.20) 

-0.10 

(1.27, 

-154.50) 

0.90 

(0.93, 

-172.66) 

-0.05 

(1.20, 

-154.77) 

0.95 

(0.93, 

-173.07) 

1.00 

(0.94, 

-173.44) 
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TABLE  B.2  (Continued) 


Values 

of  G+(kx) 

for  kb  *  5.0 

X 

C+(kx) 

X 

Out) 

-1.00 

(1.96, 

-170.85) 

0.00 

(0.75, 

143.24) 

-0.95 

(0.02, 

-30.95) 

0.05 

(0.81, 

146.42) 

-0.90 

(1.12, 

134.96) 

0.10 

(0.85, 

149.28) 

-0.85 

(1.72. 

161.53) 

0.15 

(0.38, 

151.35) 

-0.30 

(1.95. 

-176.90) 

0.20 

(0.91, 

154.15) 

-0.75 

(1.94, 

-158.66) 

0.25 

(0.94, 

156.21) 

-0.70 

(1.79. 

-142.84) 

0.30 

(0.95, 

158.06) 

-0.65 

(1.57. 

-128.94) 

0.35 

(0.97, 

159.71) 

-0.60 

(1.30. 

-116.59) 

0.40 

(0.98, 

161.20) 

-0.55 

(1.02, 

-105.55) 

0.45 

(0.99, 

162.53) 

-0.30 

(0.75. 

-95.63) 

0.50 

(1.00, 

163.73) 

-0.45 

(0.49. 

-36.70) 

0.55 

(1.00, 

164.80) 

-0.40 

(0.26. 

-73.64) 

0.60 

(1.01. 

165.77) 

-0.35 

(0.06, 

-71.36) 

0.55 

(1.01, 

166.64) 

-0.30 

(0.12. 

115.22) 

0.70 

(1.01. 

167.43) 

-0.35 

(0.27, 

121.17) 

0.75 

(1.02, 

163.14) 

-0.30 

(0.41, 

126.54) 

0.80 

(1.02, 

168.79) 

-0.15 

(0.52, 

131.39) 

0.35 

(1.02, 

169.38) 

-0.10 

(0.61, 

135.76) 

0.90 

(1.02, 

169.92) 

-0.05 

(0.69, 

139.70) 

3.35 

(1.02, 

170.41) 

1.00 

(1.02, 

170.35) 

Values 

of  G+(kac) 

for  kb  »  6.0 

X 

G+(kx) 

X 

(tat) 

-1.00 

(1.34, 

178.59) 

0.00 

(1.39, 

171.89) 

-0.95 

(0.54, 

105.91) 

0.05 

(1.36, 

173.94) 

-0.90 

(1.57, 

148.39) 

0.10 

(1.33, 

175.62) 

-0.35 

(1.31, 

179.63) 

0.15 

(1.31, 

176.98) 

-0.30 

(1.61, 

-155.21) 

0.20 

(1.23, 

178.07) 

-0.75 

(1.21, 

-134.08) 

0.25 

(1.26, 

178.93) 

-0.70 

(0.74, 

-115.94) 

•>.30 

(1.24, 

179.61) 

-0.65 

(0.27, 

-100.15) 

0.35 

(1.22, 

-179.87) 

-0.60 

(0.15, 

93.70) 

0.40 

(1.20. 

-179.-7) 

-0.55 

(0.51, 

105.90) 

0.45 

(1.19. 

-179.17) 

-0.50 

(0.79, 

116.67) 

0.50 

(1.17, 

-173.95) 

-0.45 

(1.02, 

126.17) 

0.55 

(1.16. 

-178.79) 

-0.40 

(1.13, 

134.54) 

0.60 

(1.15, 

-178.68) 

-0.35 

(1.29, 

141.90) 

0.65 

(1.14, 

-178.60) 

-0.30 

(1.37, 

148.33) 

O.'O 

(1.13, 

-173.56) 

-0.25 

(1.41, 

153.92) 

0.75 

(1.12, 

-178.53) 

-0.20 

(1.43. 

158.76) 

0.30 

(1.11. 

-178.53) 

-0.15 

(1.43, 

162.90) 

0.35 

(1.11. 

-178.53) 

-0.10 

(1.43, 

166.43) 

0.90 

(1.10, 

-178.54) 

-0.05 

(1.41, 

168.40) 

0.95 

(1.09, 

-178.56) 

1.00 

(1.09, 

-173.59) 
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TABLE  B.2  (Continued) 


Values 

of  G^(lcx) 

for  kb  -  6.2832 

X 

G+(kx) 

X 

K 

Out) 

-1.00 

(1.38, 

176.73) 

o.oo 

(1.41, 

-180.00) 

-0.95 

(0.71, 

109.01) 

0.05 

(1.37, 

-178.46) 

-0.90 

(1.71, 

153.38) 

0.10 

(1.33, 

-177.30) 

-0.85 

(1.33, 

-174.04) 

0.15 

(1.30, 

-176.46) 

-0.30 

(1.49. 

-147.33) 

0.20 

(1.27, 

-175.87) 

-0.75 

(0.97, 

-125.87) 

0.25 

(1.24, 

-175.48) 

-0.70 

(0.41, 

-107.05) 

0.30 

(1.21. 

-1,5.26) 

-0.65 

(0.11, 

39.23) 

0.35 

(1.19, 

-175.15) 

-0.60 

(0.55, 

103.55) 

0.40 

(1.17, 

-175.14) 

-0.55 

(0.90, 

116.07) 

0.45 

(1.15, 

-175.19) 

-0.50 

(1.17. 

127.05) 

0.50 

(1.14, 

-175.23) 

-0.45 

(1.36. 

136.68) 

0.55 

(1.13. 

-175.41) 

-0.40 

(1.48, 

145.08) 

0.60 

(1.12, 

-175.55) 

-0.35 

(1.55, 

152.38) 

0.65 

(1.11. 

-175.70) 

-0.30 

(1.58, 

158.63) 

0. 70  . 

(1.10, 

-175.36) 

-0.25 

(1.58, 

164.05) 

0.75 

(1.09, 

-176.01) 

-0.20 

(1.57, 

168.59) 

0.30 

(1.08, 

-178.17) 

-0.15 

(1.54, 

172.38) 

0.35 

(1.08, 

-176.32) 

-0.10 

(1.50. 

175.50) 

0.90 

(1.07, 

-175.46) 

-0.05 

(1.46, 

178.01) 

0.95 

<1.07, 

-176.60) 

1.00 

(1.07, 

-176.73) 

Values  of  C^Cot)  for  kb  ”  T.O 


X 

(Vac) 

X 

K 

(kx) 

-1.00 

(2.06, 

174.39) 

0.00 

(1.23, 

-159.46) 

-0.95 

(1.19, 

119.85) 

0.05 

(1.17, 

-160.06) 

-0.90 

(2.05, 

169.13) 

0.10 

(1.12, 

-161.01) 

-0.35 

(1.76, 

-154.76) 

0.15 

(1.08, 

-162.18) 

-0.30 

(1.01, 

-125.77) 

0.20 

(1.05, 

-163.45) 

-0.75 

(0.17, 

-101.56) 

0.25 

(1.02, 

-164.74) 

-0.70 

(0.58, 

99.10) 

0.30 

(1.01, 

-165.98) 

-0.65 

(1.13, 

116.93) 

0.35 

(0.99, 

-167.15) 

-0.60 

(1.60, 

132.39) 

0.40 

(0.98, 

-168.23) 

-0.55 

(1.37, 

145.30) 

0.45 

(0.98, 

-169.20) 

-0.50 

(2.00, 

157.41) 

0.50 

(0.97, 

-170.07) 

-0.45 

(2.04, 

167.37) 

0.55 

(0.97, 

-170.34) 

-0.40 

(2.01, 

175.82) 

0.60 

(0.97, 

-171.53) 

-0.35 

(1.94, 

-177.14) 

0.65 

(0.97, 

-172.14) 

-0.30 

(1.34, 

-171.42) 

0.70 

(0.97, 

-172.68) 

-0.25 

(1.72, 

-166.93) 

0.75 

(0.97, 

-173.16) 

-0.20 

(1.61, 

-163.58) 

0.80 

(0.97, 

-1*3.58) 

-0.15 

(1.49, 

-161.29) 

0.35 

(0.97, 

-173.97) 

-0.10 

(1.39, 

-159.93) 

0.90 

(0.97, 

-174.31) 

-0.05 

(1.30, 

-159.37) 

0.95 

(0.97, 

-1*4.61) 

1.00 

(0.97, 

-174.39) 

